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ABSTRACT 


This study presents the basic calculations from which the 
feasibility of a scanning microwave radiometer system for mapping 
the thickness of lake ice may be evaluated. An analytical model for 
the apparent brightness temperature as a function of ice thickness 
has been developed, and elaborated to include such variables as 
galactic and atmospheric noise, aspect angle, polarization, temperature 
gradient in the ice, the presence of transition layers such as snow, 
slush, and water, increased loss due to air inclusions in the ice 
layer, and the presence of multiple ice thicknesses within the antenna 
footprint. It was found that brightness temperature measurements at 
six or seven frequencies in the range of 0.4 to 0,7 GHz were required 
to obtain unambiguous thickness estimates. A number of data processing 
methods were examined; a minimum-distance algorithm and a ternary 
brightness temperature quantization scheme were found to give accuracies 
of about 1 cm for ice thicknesses up to one meter. A brief system 
study of the effects of antenna beamwidth, scanning rate, receiver 
bandwidth, noise figure, and integration time showed that the receiver 
sensitivity requirements could be met with state-of-the art components. 
The radiometric system was found to fail for sea ice thickness measure- 
ment due to the very high losses of sea ice. 
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CHAPTER I 


INTRODUCTION 

The ability to map the thickness of lake ice from an aircraft 
on a routine basis would have important scientific and economic 
consequences. For example, the shipping season on the Great Lakes 
could be extended if detailed ice thickness information were 
available on a daily basis. For a number of reasons, some of which 
will be discussed in the following paragraphs, the microwave radio- 
meter, Peake (1968), appears to be the appropriate remote sensor for 
this task. The purpose of this report is to estimate the performance 
of a microwave radiometer system over a wide range of frequencies, 
and realistic ice conditions with a view to establishing the 
feasibility and design parameters .of an operational instrument for 
the remote sensing of lake ice thickness. 

The particular choice of microwave radiometer is based on a 
number of constraints imposed by the problem itself. In the first 
place, any useful system must be independent of the weather over 
the ice surface, thus precluding the use of sensors operating in 
the optical or infrared portions of the spectrum. Secondly, an 
imaging (i.e., scanning) sensor must estimate ice thickness of the 
order of a meter or less in the presence of surface roughness 
which may be a significant fraction of thickness. This requirement 
precludes the use of the typical active microwave sensors such as 
the short pulse on carrier, the video pulse, and the side-looking 
airborne radar (SLAR). 

As for the al timeter-type short-pulse radars, these have 
proved to be quite successful in determining the properties of 
glacier or polar ice with thicknesses of the order of kilometers, 
Harrison (1973), The technical problem of producing a pulse modulated 
carrier with the pulses a few nanoseconds in duration and detecting 
the reflection from the first surface (air - ice) is not difficult. 
Such radars have been used successfully to study the wave 
structure of the ocean, Yaplee (1971). The problem in adapting to 
lake ice measurements is that of detecting the second surface (ice - 
water) in the presence of the delayed scattering from the rough first 
surface. In addition to requiring a smooth first surface, 
detection requires that the second surface be fairly well defined. 

For this case (i.e., well defined first and second surface layers), a 
successful profiling thickness measuring system has been designed 
and tested. Cooper (l976), Vickers (1974), using an S-band carrier 
frequency. A more sophisticated carrier system at 3 GHz has been 
described, lizuka (1971), with an antenna focussing system. An 
FM radar for ice and snow measurement has been described by Page 
(1973). However, a gradual transition from ice to water via a mixture 
layer of slush would act as an impedance match over several wave- 
lengths of the carrier, even though the thickness of the transition 
layer was shorter than the pulse length. 
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A more attractive short-pulse technique is that in which the 
significant spectral content of the pulse lies in the region for 
which the wavelength is comparable to the ice thickness, i.e., 
the video pulse system. Early versions, Meyer (1966), of such 
systems appeared to give ambiguous results even over smooth ice. 
However, modern versions, Moffatt (1976), Campbell (1974), with 
sophisticated antenna designs and data processing procedures should 
provide a capability for direct measurements of ice thickness. 

Because these active systems must operate at normal incidence, they 
provide a profile of the thickness along the flight path, but 
cannot scan to provide an image. A comparative study of a number 
of active systems is given by Chang (1976). 

Finally, the side-looking airborne radar (SLAR) , while it 
can provide an excellent all weather image of the extent and 
structure of the ice cover, does not seem adaptable to quantitative 
measurements of ice thickness. Even if the gray scale of the image 
was calibrated quantitatively, the observed intensity, particularly 
for the near grazing angles at the edge of the image, is determined 
primarily by the first surface scatter. Moderately successful 
attempts have been made to distinguish various types of ice, but 
the thickness per se could not be estimated at all. Rouse (1969), 
Larrowe (1971), Parashar (1974). Thus, active microwave sensors, 
for a variety of reasons, are not appropriate for producing images 
of the distribution of ice thickness. 

The passive microwave (radiometer) systems, on the other hand, 
because their primary response is controlled more by dielectric 
constant than by roughness, appear to offer a good possibility 
for ice thickness measurement. In a qualitative way the dramatic 
contrast between the brightness temperature of sea ice (with a 
dielectric constant of three and an emissivity of over 90%) and 
water (with emissivities typically of 35% in the microwave region) 
may be seen in the images reported by, Lade>y (1967) at X-band, 
and by Catoe (1967) at 19 GHz and in the UHF band by Adey (1973). 

More recently, the routine processing of data from 19 and 35 GHz 
images in a false color presentation have resulted in the ability 
to identify and map several categories of sea ice (but not ice 
thickness) on a synoptic basis, Campbell (1974), Gloersen (1974). 
Since the quantitative thickness measurement depends on having the 
electrical wavelength comparable to the layer thickness, these 
higher microwave frequencies are clearly more suited to layers a 
few millimeters thick, as one may see from the successful profiling 
of oil spill thicknesses reported by Hoi linger (1973), and by 
laboratory measurements of ice layers a few millimeters thick 
reported by Edgerton (1971). More recently, several preliminary 
attempts have been made to demonstrate the ability of the radiometer 
as a sea ice thickness sensor in the profiler mode, Adey (1972), 

Tiuri (1975), using significantly lower frequencies (e.g., 600 MHz). 

A number of other researchers have experimentally observed the inter- 
ference phenomenon of layered media, England, Pascalar (1964), 
Basharinov (1971), Jean (1971), Blinn (1972), Schmugge (1973), 
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Popov (1976). Thus, measurement of layer thicknesses by radiometric 
methods may be considered established, but it is clear that the 
design of a system capable of producing an image of the spatial 
distribution of ice thickness requires a thorough study before a 
suitable system configuration can be chosen and the design 
parameters established. 

This report presents the basic calculations from which the 
feasibility of such a system may be evaluated. A number of 
figures are used to exhibit the apparent brightness temperature 
of layers of ice as a function of ice thickness, frequency, 
polarization, angle of incidence, surface roughness, etc.. Several 
different multi -frequency measurement systems and data processing 
methods are also considered. However, before presenting the 
details of these parametric calculations, it may be useful to 
provide a qualitative discussion to place the thickness measure- 
ment problem in perspective. (For similar discussions from a 
somewhat different point of view, see England, Popov (1976)). 
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CHAPTER II 


GENERAL FEATURES OF A LAYERED MEDIUM RADIOMETRIC RESPONSE 

The brightness temperature T, of a plane layered medium at a uniform 
physical temperature T^- (see Figure 2) is composed of two terms 




(2-1) 


of which the first, Tu , represents the brightness temperature of the 
radiation emitted by the layers, and the second, T^^, represents external 
radiation reflected by the layers. These two terms may be found from the 
emission coefficient e of the surface. 


Tbe = ® T,. 


Tbr= n-e)(Tq + V 


(2-2a) 


where (1-e) = r is the power reflection coefficient of the surface, T 
is the brightness temperature of the incident cosmic radiation and T^^ 
is the brightness temperature of the radiation emitted by the atmosphere. 


To illustrate the general features of the response of a 
radiometer viewing a layer of ice over an (effectively) infinite 
depth of water, it is sufficient to know the voltage reflection 
coefficient R for an incident plane wave. 


R 


- 2vd 

- R 2 e ^ 

1 - R^ R 2 e’ 


(2-3) 


Rl = (1 - /ip/(l + /T^-) 

R2 = (•% - ^p/( + /ip 


where 

R-j = voltage reflection coefficient (air to ice):ii-2/7 

R 2 = voltage reflection coefficient (water to ice)lL+2/3 

d = thickness of ice 

Y = a + j3 = complex propagation constant in ice 
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A^. = wavelength in ice 

3 = 2tt/a^ propagation constant in ice 

a = attenuation constant in ice 
e^. = complex relative dielectric constant of ice 

= complex relative dielectric constant of water 


For angles of incidence other than normal, one need only replace 
R-i > Rp Fresnel reflection coefficients appropriate to the 

incident polarization state, and replace d by ii, the path length 
in the ice, i = d/(l - sin^ Note that the rate of 

oscillation of the reflection°coefficient, given by the phase 
23^, is independent of polarization and only weakly dependent on 
angle of incidence. The magnitude of the oscillation is controlled 
by the reflection coefficients, and these are also almost independent 
of look angle for favorable polarization states (circular polarization 
or 45° linear). 

If, for the purposes of a qualitative discussion, one assumes 
that the reflection coefficients R*j, R 2 are real, then the emission 
coefficient of the surface, e = 1 - r has the form 


_ (1 - r^-) (1 - A r^) 

(1 + A r. r. + 2 JK r. r, cos 23 d) 

I W I w 


(2-4) 


*^i ~ ^i ^i power reflection coefficient (air to ice)2.2/25 
r^ = R^ R^ power reflection coefficient (ice to water)_;^ 4/9 
r = 1-e = RR* power reflection coefficient of layers. 

A = exp(- 4a d). 


To illustrate the behavior of the reflection coefficient. Figure 1 
shows T^e as a function of d/A^ (the ice thickness in wavelengths) for 
various values of the parameter a/3. (An attenuation of 1 dB/wavelength 
in ice corresponds to a/3 = (407 t log^jQ e)"* * .0183.) 

For fresh water ice, with attenuation of the order of 1 to 
3 dB per meter, (approximately independent of frequency) in the 
lower microwave range, the curve a/3 = 0 illustrates the general 
behavior of the brightness temperature dependence. For sea ice, 
with much higher attenuations, the curves a/3 = 0.1 or 0.3 may 
be more representative. 
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Note that these curves represent ice with constant attenuation 
per wavelength (i.e., a proportional to frequency), whereas both 
sea ice and fresh water ice tend to have a independent of frequency. 
It is significant, however that even for the case a /3 = 0.3, repre- 
senting an attenuation higher than any likely to be encountered in 
fresh water ice, the period of the oscillations of with thickness 
is essentially unchanged. 

The nature of these oscillations may be placed in better per- 
spective by writing Equation (2-4) as a Fourier cosine series 


T.. (1 - r ) (1 - A r ) 
T = -J ] ! 

0 - A r r ) 


[1 - 2 (r. cos 2 N 


( 2-5 a) 

+ 2 (r. A) cos 4g d -2(r. cos 6gd + ... 

= Z cos 2n8d n = 0,1,2 ... (2-5b) 

which, for e-j = 3.1, = 78, a = o becomes, as an illustration 


T|:^= 145 - 53 cos 23d + 10 cos 43d - ••• °K. (2-6) 


Similarly, the reflected contribution is given by the corresponding 
series 


■’’hr = ( VTa) (1 - ^ '"w) [A + 2(r. r cos 26 d 

"" 3 (1 - A n r„) ° ' 

- 2{r- r^ A) cos 43 d + . . . ] 

= I T^' cos 2n3d n = 0,1,2 ... (2-7) 

*0 " (>"i ^ "“w ■ ^ ^ ''i ^ 

which, for Tg + Ta "" 1000°K becomes, as an illustration, 

Ttjy, = 470 + 200 cos 23 d - 37 cos 43 d + . , . (2-8) 
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In the series form, the effect of changes in attenuation, finite 
instrument bandwidth, and fluctuation in ice thickness are easily 
found. For example, if the radiometer has an ideal band pass filter 
characteristic, passing frequencies corresponding to the range +^3/2, 
and the ice thickness within the antenna footprint is uniformly 
distributed in the range d ± Ad/2, then the averaged Tj^g becomes 
approximately, for ice of fairly low attenuation. 



00 



\ 

1 \ 

sin ndA3 

sin n3Ad 

ndAg 1 

, n3Ad 

/ 

^ / 


cos 2n 3d 


(2-9) 


i.e., the amplitude of each harmonic is simply reduced by the two 
"sin x/x" factors; larger attenuations will produce further small 
changes in the coefficients over and above the explicit factor 
exhibited in Equation (2-3a). Similar effects due to random 
fluctuations in dielectric constant are discussed in Tsang (1976). 

A more detailed discussion of the effects of finite bandwidth for 
various bandpass shapes, is given in Popov (1976). The effects of 
a finite beamwidth may also be found by averaging over the range of 
path lengths ^( 0 ) determined by the range of look angles in the beam. 


The behavior of the total perceived brightness temperature Tj^ 
will depend on the relative magnitudes of the temperatures (Tg+Tg) 
and T. 270°K. At low frequencies, below about 150 MHz, the^cosmic 
noiseHemperature Tg is much larger than Tj , and strongly frequency 
dependent Ko, (1958). Typically, the major portion of the radio 
sky has temperatures larger than Tg > 2/f2*' (f in GHz); for fre- 
quencies around 0.1 GHz this may conveniently by written 
Tg ^ 1000/(f/* 1 )2‘^. Thus below 150 MHz or so T^ is dominated 
by the reflected term 


Tb Ai T^y. " I Tg 


cos 2 and 


Tn' 


from (2-7) 


and since the signs of Tg' are opposite those of Tg (n<0), the 
oscillation of T^y. will be 180° out of phase with those shown in 
Figure 1 . 

Above 350MHz on the other hand, Tg < 80° K over most of the radio 
sky and the brightness is dominated by emission. 


Tfa i Tjjg = ^ T„ cos 2 end . from (2-5) 

At intermediate frequencies, between 150 MHz and 350 MHz, as will be 
seen subsequently (see, e.g., Figure 6) the two components of bright- 
ness are comparable in magnitude, and the variation of brightness with 
ice thickness depends-'on the (usually unknown) ratio Tg/T^-. 
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From the perspective of the two cosine series (2-5b) and (2-7) 
for Tbr and T\^q, the approach to a thickness measuring system becomes 
evident. At one extreme, (see Figure 1) it is possible to operate 
in the first half cycle of cos 2 ad, thus avoiding any problems of 
ambiguity by restricting the phase to lie in the range 0 < 2 ad < tt. 
For ice of one meter thickness, this requires f < c/(4 d/^y) -fb 42 MHz. 
At these frequencies, the brightness T|^ is completely dominated by 
the reflected cosmic noise. Since the cosmic noise will vary (it is 
dependent on the direction of the flight path and angle of incidence) 
one must either measure Tg with an upward looking antenna, or avoid 
the problem by making certain assumptions about the behavior of Tg. 

For example, if it is assumed that Tg has a power dependence on fre- 
quency, Tg = c/f”^ where c is a constant, one might estimate the 
thickness by observing the reflected brightness temperature and its 
derivative 



dTbr 


\ 


I f I \ df j 


-(Tq'/T,-' + Cos j.) 

m(TQVT-* + Cos i|j) + Sin 


To’/Ti' =A^/2/?ri^ 

ip = 2$Ji (2-10) 


This transcendental equation is now independent of Tg, and will 
have unambiguous solutions if is somewhat less than 7T,^i.e., the 
system would be restricted to frequencies below about 42 MHz. More 
eleborate systems, operating at somewhat higher frequencies ('^ 100 MHz) 
are discussed in Chapter III. These systems reduce, but do not 
eliminate the effects of ambiguity by utilizing two or three frequencies. 

The more straightforward method for estimating ice thickness is 
made evident by Equation (2-5) and Figure 1, both of which indicate 
that the period of oscillation of T^ ( 23 d) is changed only slightly 
by attenuation, instrument bandwidth® or small variations in ice 
thickness, although both the "dc" level of the oscillations and 
their amplitude may vary. To extract the thickness, it is only 
necessary to sample the brightness at a sufficient number of 
different frequencies to establish the period of 23^. The number 
of frequencies required, and their separation, may be estimated 
from the sampling theorem, and the required maximum and minimum 
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ice thickness of interest and respectively. 

To establish the frequency interval between channels, note 
that the sampling theorem requires that Tu(23d) be sampled 
twice per period for the thickest ice to be measured. Thus, the 
frequency interval df between channels must be 

df < c/(/I7 4 c = 3 X 10^ m/sec (2-11) 

I iflaA 

i.e., less than 43 MHz for ice with dj^g^^ = 1 meter. On the other 
hand, because at least one half period must be established for 
the thinnest ice, the total span of frequencies Df must be 
greater than 


Df > c/(/S7 4 d^^. J (2-12) 

and the number of frequencies required is N ^ ^'^max^'^in * 

The choice of the highest usable frequency will be determined 
by the uniformity in ice thickness, because the amplitude of oscilla 
tion of TL(2gd) is diminished by the factor sin (3Ad)/(3Ad). 

Thus, if one assumes that 3Ad = 2 tt/3 = x (so that sin x/x = .4), 
then the highest usable frequency f^^gx given by 


^max 


(2-13) 


Thus, if A(i is, say, 10 cm, then fp^gx "" lower limit, 

in practice, will be determined by problems of antenna size, and 
by the fact that below about 350 MHz, the magnitude and variability 
of the reflected cosmic noise significantly influences the perceived 
brightness temperature of the ice layer. 

Thus, the basic "brute- force" design for an ice thickness 
measuring system would consist of a circular or 45° linear polarized 
multi frequency radiometer with N ^ (1 + W/^min^ channels, with 

a frequency separation between channels of df = c/(4 cimax) ^ 
highest frequency less than ffpgx < c/(3 Ad), and a lowest 

frequency that should be somewhat above > 300 MHz. For an 
airborne imaging system, the image can be generated by scanning 
across the flight path, or by utilizing simultaneous multiple 
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beams. In the former case (i.e., scanning, which presents the 
more severe requirement on system sensitivity), it is not hard to 
show that if the angular resolution (antenna beamwidth) is 0 ^ 
radians and the system sensitivity (precision in measuring bright- 
ness temperature) is Alg °K, then a system figure of merit is given 
by the product, Peake (T968), 


0 


B 



'in 4^ 


( 2 - 14 ) 


where 

B 

V 

H 


F = receiver noise figure with = 290° K 
= antenna temperature -iMg 

For an aircraft with V/H*^.08 sec“\ B = 5 MHz, - 500° K, 
for example, 

ATg = 1/8 (°K-radians) . 

Such a radiometer with, say, 10° angular resolution (en'^l/O) 
would allow a precision in measuring Tn of aTd'^.75° k , while a 
system with 3° angular resolution (0p.-::rl/2O) would have 
AT^ = 2.5° K. 

It will appear subsequently that errors in measuring Tn of 
as much as 5° K will not seriously affect the ability of a System 
to estimate ice thickness. From the superficial view of system 
design implied by Equation (2-8), it is clear that the limited 
dwell time imposed by the scanning requirement would not become 
a problem (to be solved, for example, by multiple simultaneous 
beams) except for systems with high resolution (beamwidths less 
than a few degrees) or with large V/H ratios. Furthermore, the 
required radiometer system performance is well within the current 
state-of-the-art for radiometer hardware. Thus, the major technical 
problem in producing an imaging, thickness measuring radiometer is 
the development of a flyable, low loss, low VSWR, multichannel, 
scannable (or multi beam) antenna in the 300 MHz to 800 MHz frequency 
range. The antenna design problem is not considered in this report, 
the remainder of which is devoted to justifying the statements made 
in this brief survey, and to developing a system which will produce 
valid ice thickness estimates with as small a number of radiometer 
channels as possible. 


= bandwidth of individual radiometer channel (Hz) 
= vehicle velocity (m/s) 

= vehicle altitude (m) 

= (F - 1) T + T„ 
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CHAPTER III 


APPARENT BRIGHTNESS MODEL OF A LAYERED MEDIUM 
0. Introduction 


In this section, the theory for calculating the apparent 
brightness temperature of a layered medium is reviewed and used 
to calculate the expected brightness of a smooth ice layer over 
water. In order to introduce the several physical parameters 
controlling the perceived brightness, the case of normal incidence 
with ice and water at the same temperature is considered first. 

In the following chapter, the model is elaborated to include the 
effects of transition layers of various types. 

1 . Brightness Model Geometry 

The basic problem geometry is shown in Figure 2. An infinitely 
deep layer of fresh water is assumed to be covered by a layer of 
ice of uniform thickness d. The ice temperature T. and water 
temperature T are both assumed to be 273° K. The^ surface layer 
is illuminated by galactic noise T and atmospheric noise T . 

The observing antenna will then receive two contributions from 
the direction of the surface, namely the reflected sky noise temper- 
ature and the noise temperature emitted by the ice and water. The 
total apparent brightness temperature T. of the radiation incident 
on the antenna is then given by ° 

= (1 - |Rl^) T. + |R|^ (T + Tg) (3-1) 


where emission (Tj^^) reflection (T^^,) 

R is the voltage reflection coefficient of the ice-water 
combination 

T^. is the physical temperature of the ice and water 
Tg is the brightness temperature of the galactic noise 
T^ is the brightness temperature of the atmospheric noise. 
For normal incidence, 

- z 

R = — 2. (3_2) 

Z, + Z 


where 


12 




Figure 2. The apparent brightness temperature of the ice layer/ 
water system a is the sum of the emitted component 
Tbe=(l - |Rr) T^. , and the reflected component T|^^= 

|R|2(Tg+T^). 
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* 377 

Z + Z. tanh yd 
w 1 _ 1 

1 . + 2 t^aTili^Y^d 
1 w M 




Y . = 

’ 1 C 1 


and 


is the relative complex dielectric constant of ice 
is the relative complex dielectric constant of water 
d is the thickness of the ice layer 
Z^ is the free space impedance 
Z^ is the ice impedance 

Z^ is the water impedance 

is the complex propagation constant in ice 
c is the free space propagation velocity 
This basic model assumes: 

- the apparent brightness temperature is observed at normal 
incidence 

- the ice layer has a smooth surface on the top and bottom 

- the top surface of the ice is free of snow, slush, or water 

- the bottom surface of the ice has no slush at the ice-water 
interface 

- the ice layer has no air bubbles in it 

- the water and ice are at the same uniform temperature. 

The effect of relaxing each of these model constraints will be 
considered in the following chapter. Before doing so, the input 
parameters r . e . , T^, and T, will be defined. 

wig d 
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2. Brightness Model Parameters 
A. Water Model 

The dielectric constant of water may be represented by an 
equation of the Debye form. The particular model used in this 
report is the one developed by Stogryn (1971). 

4 . ^dw . g 

% ^ ^«w 1 + j 2 ITT f ^ t V f (3-3) 


where 


e is the high frequency dielectric constant 

<xyj 

S is the salinity in parts per thousand 
T is the temperature in °C 
N is the normality of the solution 
is the static dielectric constant 
g is the ionic conductivity of the dissolved salts 
T is the relaxation time in seconds 
f is the frequency in Hz 


N = S(1.707 X 10"^ + 1.205 x 10"^ S + 4.058 x 10‘® S^) .9141 

a = 1.0 - 0.2551 N + 5.151 x 10"^ - 6.889 x 10'^ 

e, = a(87.74 - 0.4008 T + 9.398 x 10'^ + 1.410 x 10‘® T^) 

A = 25.0 - T 

g = 2.003 X 10'^ + 1.266 X 10‘^A + 2.464 x 10‘® i? - 

5(1.849 X 10"® - 2.551 x lO'^A + 2.551 x 10"® A^) 
o = 5(0.182521 - 1.46192 x 10"^ S + 2.09324 x 10'® 5^ - 
1.28205 X 10"^ 5^) exp (-Ag ) 
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b = 0.1463 X 10"^ N T + 1.0 - 0.04896 N - 0.2967 + 

5.644 X 10'^ 

2ht = b(1.1109 X 10"'° - 3.824 x 10'^^ T + 6.938 x 10"^^ - 

5.096 X 10"’^ T^) 

= 8.854 X 10 farads/m 

These equations are for sea water of standard composition. A recent 
study by Klein and Swift shows that for a NaCl solution, the factor 
,9141 in the expression for the normality N should be replaced by 
uni ty. 

B. Ice Model 

The model for the dielectric constant of the ice layer, as 
for the water, is of the Debye form. The main difference between 
the two is that the relaxation frequency for ice is in the KHz 
range and for water in the GHz range. The relative complex 
dielectric constant of fresh water ice is given by 


e 


i 



^di “ ^°°i 


where 


is the high frequency dielectric constant 
is the static dielectric constant 
f^ is the relaxation frequency in Hz 
T is the temperature in °K 


e^.j = 90. - (T - 273.) 0.3581 
e . = 2.846 + 0.001333 T 
f = a« T^ + a -3 T^ + a« T T 

a. = 0.11666643 x 10"^ 
a 3 = -0.11573310 x 10*^ 
a 2 = 0.43053546 x 10^ 
a^ =-0.71170619 x 10^ 
a = 0.44104997 x 10® 


(3-4) 
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These parameters were derived by fitting the Debye model to 
measurements of Hasted (1961), and Iglesias (1967). A comparison 
of the loss tangent given by this ice model with the measurements 
upon which it is based is shown in Figure 3. The real part of the 
relative complex dielectric constant shows no dependence on 
frequency but does decrease from 3.21 to 3.18 as the temperature 
decreases from 0°C to -40°C. This ice model corresponds to very 
good ice, i.e., ice without air bubbles, cracks, etc.. At 
frequencies well above the relaxation frequency (a few KHz)^ this 
model produces a loss tangent that is approximately proportional 
to wavelength so that the attenuation is roughly frequency 
independent, (a ^ 163 * tan 6 * f in dB/meter with f in GHz.) 

The values of attenuation recently measured for actual lake ice^ 
Vickers (1975), ranged from 1 to 3 dB/meter (depending on the 
condition of the ice), almost independent of frequency in the range 
1 to 10 GHz; these attenuations correspond to 

.006/f < tan <.019/f with f in GHz. 

C. Galactic Noise Model 

The noise temperature contribution T from the galaxy is 
variable with time, frequency, and direction. For the purposes 
of the calculations made here, it will be assumed that the cosmic 
background is isotropic but of variable intensity with a (-2.7) 
power dependence on frequency. The model used is 

T = G/f^-^ (3-5) 

9 


where 


G is the galactic temperature factor (2 G < 40) 
f is the frequency in GHz. 

The bounds of possible galactic temperatures are shown in Figure 4. 
The actual temperature distribution in the sky has a complex 
structure, Ko (1958), including individual radio sources. Figure 5 
shows typical radio maps of the sky at several different frequencies. 
The effective T is thus not only a function of position in the sky 
but of the anteSna beamwidth as well. For this reason^ Figure 4 can 
be considered only as an approximation for the temperature which 
might be encountered. The maximum temperature (G = 40) would be 
observed in the direction of the galactic center, but the brightness 
falls off rapidly toward the minimum (G = 2) as the look direction 
moves toward the direction of the galactic pole. The solid angle 
subtended by any given brightness temperature contour also decreases 
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Figure 3. The loss tangent of an ice model fitted to measurements presented by 
Iglesias. 
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Figure 4. The extremes of galactic noise temperature which may 
be encountered. 








rapidly with frequency. Thus, to evaluate precisely the cosmic 
noise contribution for any particular mission, it would be necessary 
to take into account the actual antenna pattern and the flight 
path, relative to the galactic coordinates. 

A general survey of the cosmic noise distribution over the 
entire range of frequencies of interest here is given by Ko (1958). 

The sky maps of Figure 5 are taken from this reference. More 
detailed brightness maps are given by Seeger (1959) and Landecker 
(1970) at 400 MHz and 150 MHz respectively. 

D. Atmospheric Noise Model 

The atmospheric noise model is based on data published by Hogg 
(1961). The atmospheric contribution to the sky temperature at low 
frequencies (1 GHz or less) shows little dependence on humidity, 
and consequently, all calculations are made assuming the standard 
atmosphere. In the frequency range 0.1 GHz to 1.0 GHz, the 
atmospheric noise contribution is on the order of 5°K or less. 

E. Urban Noise 

In addition to the universally present cosmic and atmospheric 
noise, there are locally severe sources of man-made noise which 
may influence radiometer operation. The coherent sources of man- 
made radiation include radar, television, communications equipment, 
navigation aids, etc.. These sources of interference are a function 
of frequency band and geographic location and will not be considered 
in detail here. The incoherent noise is due to such sources as 
ignition systems, power lines, and electric motors. Investigation 
of this type of noise has been reported by a number of authors, 

Ploussios (1968), Skomal (1969). It has been found that the 
measured temperatures over large urban areas may be quite high, 
many thousands of degrees at 369 MHz for example, but that the 
brightness decreases rapidly beyond the limits of urbanization. 

For example, the brightness temperature dropped from 48000 °K to 
0 / 300 ®K in one mile as an aircraft flying at 2000 feet passed 
over Newark, New Jersey, Ploussios (1968). Thus, an aircraft flying 
at a few thousand feet, a mile or more off-shore would not be 
expected to experience interference from incoherent urban noise. 

A special investigation, however, would have to be made to ensure 
that the chosen radiometer frequencies do not coincide with possible 
sources of coherent interference. 

3. Brightness Model Power Reflection 

Coefficient and Temperature 

As can be seen from Equations (3-1) and (3-2), the apparent 
brightness temperature depends on the reflection coefficient, and 
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this in turn is a function of the ice layer thickness d. A 
representative plot of the reflection coefficient as a function 
of d for two different frequencies is shown in Figure 6. The 
apparent brightness temperatures corresponding to those frequencies 
for three different levels of galactic temperature are shown in 
Figures 7(a), (b), and (c). The power reflection coefficient is 
also a sinusoidal function of frequency. Figure 8 shows the 
frequency dependence of the power reflection coefficient for two 
different ice thicknesses. The corresponding brightness temper- 
atures are shown in Figures 9 (a), (b), and (c) for several 
galactic temperatures. The model assumed here for good ice has 
very small but non-zero loss. Thus, as the thickness of the 
ice layer increases, it will tend to become more absorptive. This 
effect can be seen in Figure 8., where the oscillation in the power 
reflection coefficient for the 100 cm ice layer is somewhat smaller 
than for the 10 cm layer. Similarly, the emissivity increases with 
corresponding slight increases in apparent brightness temperature, 
as shown in Figure 9 (a), (b), and (c). 



different frequencies. 
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Figure 7 (c). The apparent brightness temperature of the ice 
layer as a function of layer thickness and 
magnitude of galactic temperature. 
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Figure 8. The power reflection coefficient as a function of 
frequency for two different ice layer thickness^. 
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Figure 9 (c). The apparent brightness temperature of the ice 
layer as a function of frequency and magnitude 
of galactic temperature for two different ice 
thicknesses. 


For low frequencies, the galactic temperatures are very 
large, the reflection mechanism predominates, and the apparent 
brightness temperature curves resemble in form those of the power 
reflection coefficient. As frequency increases, the galactic temper- 
ature decreases to the point where T + T equals T. (about 150 MHz 
for small G) , the sum of the emi tted^power and reflected power be- 
comes constant, and the apparent brightness temperature becomes 
independent of the power reflection coefficient, and hence of the 
ice thickness d. At still higher frequencies, the emission clearly 
dominates the reflection, and the character of the versus d curves 
is reversed. It can be seen from Figures 7 (a), (b) , (c) and 
Figures 9 (a), (b), and (c) that the frequency for which the 
emissivity equals the reflectivity, is a function of the galactic 
temperature T . At these "null" frequencies the apparent brightness 
temperature i§ independent of ice thickness. Clearly frequencies 
in the range 100 - 300 MHz are unsuitable for a general purpose 
instrument, because the reflected cosmic noise contribution depends 
strongly on the position of the galactic plane relative to the 
reflected beam. 
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CHAPTER IV 


APPARENT BRIGHTNESS MODEL GENERALIZATIONS 

0. Introduction 

The theory for calculating the apparent brightness temperature 
of a layered medium at normal incidence was reviewed in the 
previous chapter. Since the major advantage of the radiometer 
over active systems is the ability to scan, it is now desirable 
to include the effects of look angle and polarization. The model 
is also generalized to include the effects of temperature 
gradients in the ice. In the natural state^ there may also be 
transition layers such as snow, slush, or water on top of the ice 
layer. In addition, natural ice will have some texture on the 
top surface and may also exhibit significant thickness variations. 
Such ice will also have losses greater than those predicted by 
Equation (3-4)^ due to air inclusions, etc.. The effect of these 
conditions on the apparent brightness temperature will be cal- 
culated in this chapter. 

1 . Variation with Look Angle and Polarization 

If there existed only a single interface between two media, 
the reflection coefficient as a function of scan angle would be 
given by the well-known Fresnel reflection coefficients. When 
there are one or more homogeneous layers, the reflection and emission 
coefficients take a more complex form. Although^ closed form solu- 
tions may be easily found for the single layer case, it is more 
convenient to approach the problem by use of an iterative procedure 
described by Richmond (1965). This procedure can handle any number 
of layers, and so it will be useful for approximating the effects 
of temperature gradients in the ice. The geometry of the problem 
is shown in Figure 10. Each layer has thickness d , uniform 
thermometer temperature T , and relative complex eTectric and 
magnetic permittivities s': and u . Consider a plane wave incident 
from the direction e. For horizBntal polarization. 


= A i ® ^ 

inc h X 


(4-1) 


where A. is the incident electric field strength. In each layer 
the appropriate fields may be written 
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Figure 10. Geometry for modeling temperature gradients in the ice layer. 
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e!! = (A„ B„ e'"" 


n ' n 


(4-2) 


„ -Y (jky sins) 

K - (- ^n) ^ ^ " )« 


where 


Z = (p / e ) 

0 ^01 o' 


1/2 


(4-3) 

(4-4a) 


= jk(n E - sin^ 8)^'^^ 
^ n n ' 


(4-4b) 


s t 

By enforcing boundary conditions, the coefficients of the n + 1 
layer can be found, Richmond (1965), from those of the layer below it, 


A^i=PA+QB 
n+1 n n ^n n 


B = R A + S B 
n+1 n n n n 


(4-5a) 

(4-5b) 


where 


P„ = (.5) (1 + p^)e 
= (-5) (1 - P„) e 


= (.5) (1 - p^) e 


S„ = (.5) (1 + pj e 


' ^n+1 ^ 
n - ^n+1^ 

^ ^n+1 ^ 






(4-6a) 

(4-6b) 

(4-6c) 

(4-6d) 

(4-6e) 


The lower half space is assumed to extend to and therefore 
B =0. Setting A = 1. and applying the recursion relations 
y?elds A and B . The power reflection coefficient is 
then giv?n‘by ^ 


|Rhl^= 


(4-7) 
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For the case of vertical polarization, the results are identical 
except that now the fields are found from an assumed incident 
magnetic field 


H = A i ® 

i nc V X 


(4-8) 


with the fields in the n^*^ layer being 


= (A„ e^n"+ B„ e'^ ^ 


n n 


(4-9) 






-Y 


n ) g{jky sine) 


(4-10) 


where now 


and 

^ " l^n+l^\+l * (4-12) 

The variation of the power reflection coefficient with angle 
of incidence and polarization is shown in Figures 11 and 12. 
Vertical polarization clearly exhibits a Brewster angle 
phenomenon which causes a reduction in the dependence of the 
apparent brightness temperature on the ice layer thickness 
at a certain scan angle. These changes in the magnitude of 
the oscillation of the reflection coefficient with look 
angle, which may have an undesirable effect on certain 
thickness estimation systems, may be compensated for by 
utilizing circular or 45° linear polarization. In these cases, 
assuming that the horizontal and vertical components of the 
sky and cosmic noise are uncorrelated, the effective power 

reflection coefficient is given by ([K. + \R \^)/Z, As 

[R^I decreases with increasing 6, so increases, and the 

average of the two shows very little variation with 6. A 
plot of the power reflection coefficient for circular polar- 
ization is given in Figure 13. 


31 




0 10 20 30 40 50 60 70 80 90 100 

ICE THICKNESS (cm) 


Figure 11. The power reflection coefficient as a function of 
aspect angle and ice layer thickness for horizontal 
polarization. 
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Figure 12. The power reflection coefficient as a function of 
aspect angle and ice layer thickness for vertical 
polarization. 
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Figure 13. The power reflection coefficient as a function of 
aspect angle and ice layer thickness for 
circular polarization. 


As a check of the foregoing iterative procedure, the complex 
voltage reflection coefficient versus frequency was transformed 
to the time domain using a numerical FFT routine. (To avoid 
computational difficulties inherent in a true impulse response, 
the reflection coefficient spectrum was multiplied by the 
spectrum of very narrow raised cosine pulse.) Figure 14 shows 
the response of a 100 cm ice layer to excitation by a pulse of 
unit amplitude at normal incidence. The spacing between the 
successive pulses just corresponds to the two-way travel time 
through the ice slab and their magnitudes to the "air-ice", 
"ice-water", etc., voltage reflection coefficients. 

Although^ the motivation for this work has been to solve the 
problem of an ice layer over water, the procedure is completely 
general. Any inhomogeneous, lossy, dispersive, planar medium can 
be handled at any incidence angle or polarization. Although, 
analytic expressions for the reflection and transmission coefficient 
of certain graded dielectrics are known, the piecewise linear 
approximation has been found to be superior computationally, 
Sivaprasad (1976). Furthermore, the response to any time domain 
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Figure 14. Normalized ice layer response to narrow unit pulse 
exci tation. 


excitation can be obtained by simply multiplying the excitation 
spectrum by the complex reflection coefficient and transforming 
the product to the time domain. The response of several different 
types of media for a wide sine squared pulse excitation has been 
calculated by Simpson (1976). 


2. Temperature Gradient in Ice 

On cold days, with the air temperature well below freezing, 
a physical temperature gradient will exist in the ice layer. For 
this case, the numerical procedure given in the preceding section 
can be extended to calculate the apparent brightness temperature 
by utilizing a form of reciprocity principle. This asserts that 
if a plane wave were incident on a surface from direction e, and 
if a fraction F of the total incident power were absorbed in a 
certain layer having the physical temperature T, then the contri- 
bution of that layer to the total brightness temperature of 
radiation leaving the surface in the direction 0 would be FT. 
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The brightness temperature of a medium composed of uniform layers 
is then 



N 

I 

n=0 


F T + iRr (T + T ) 


n n 


(4-13) 


where 


= (P - 
^ n 




is the fraction of the. cower in a hypothetical incident plane 
wave absorbed in the n^" layer, and P is the net power incident 
on the layer. " 


For horizontal polarization, 


P = Re(- H^*) 
n ' n n 'z=z. 


P. = A A cose /Z 
me n+1 n+1 ' o 


For vertical polarization^ 


P = Re(E^ H^*) , 

n ^ n n 'z=z. 


P. = A A . 1 * cose Z 
me n+1 n+1 o 


(4-14) 

(4-15) 


(4-16) 

(4-17) 


To find how well a continuous, linear temperature gradient 
may be approximated by discrete, homogeneous layers, a number of 
cases were calculated. A typical result is shown in Table 1. 
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TABLE 1 

Linear Temperature Gradient Approximation 


No. 

of Layers 

f=0. 

lOGHz 


f-0.40GHz 



|R F 

T, 

be 

^b 

IR 1^ 

^be 


2 

.30927 

188.50 

500.27 

.35406 

176.29 

186.71 

4 

.30894 

188.59 

500.02 

.34729 

178.13 

188.36 

10 

.30882 

188.62 

499.93 

.34520 

178.71 

188.87 

20 

.30878 

188.63 

499.90 

.34437 

178.93 

189.07 

40 

.30876 

188.64 

499.89 

.34393 

179.05 

189.18 

80 

.30874 

188.64 

499.88 

.34370 

179.11 

189.24 

100 

.30874 

188.64 

499.87 

.34366 

179.13 

189.25 

200 

.30872 

188.65 

499.86 

.34355 

179.16 

189.27 


Total Ice Thickness = 50 cm 
Air Temperature = -40°C 
Water Temperature = 0°C 
Normal Incidence 

Galactic Temperature Factor G = 2. 

The temperature profiles were chosen such that for example, in the 
case of four layer approximation, the top 12.5 cm of the ice was 
set to -40°C, the next 12.5 cm to -30°G, the next 12.5 cm to -20°C 
and the bottom 12.5 cm was set to -10°C. 

It is clear that for fresh water ice, for which the attenuation a 
does not change much with physical temperature, only a few homogen- 
eous layers are needed to simulate a linear continuous temperature 
distribution in the ice layer. Because of the low thermal conduct- 
ivity of ice, and the fact that the temperature distribution is 
controlled by a diffusion equation, a sudden change in the air 
temperature at the surface of the ice will produce an internal 
effect only after some considerable time, proportional to the 
square of the distance from the surface. Pounder (1965). It would 
take a number of days for a linear temperature gradient to be 
established in a 50 cm ice layer for example. The actual temper- 
ature gradient in the ice would then be approximately a linear one 
reflecting the long term temperature conditions, with a daily 
variation confined to the top few centimeters of the ice layer. 

The mathematical model developed here can approximate any tempera- 
ture gradient as accurately as desired, but as may be seen from 
Table 1, changes in the physical temperature of the ice have 
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negligible effect on the apparent brightness temperature. Thus^ 
all subsequent calculations will assume a uniform ice temperature. 

3. Snow on Top of Ice 

In the natural state^it is quite possible that a snow cover 
may exist on top of the ice layer. Natural snow is very difficult 
to model accurately because of the great variation possible in the 
relative complex dielectric constant which is a function of the 
snow density and moisture content. Investigators report snow 
densities ranging in value from .06 gm/cc to .517 gm/cc, Yosida 
(1963), Kennedy (1966), Yoshino (1967), Meier (1971), Hoekstra 
(1972), Sackinger (1972), Vickers (1972). According to Gumming 
(1952), this span of densities for dry snow results in the 
permittivity ranging between 1.1 and 2.0, and the loss tangent 
(at f = 9.375 GHz, and 0°C) between .0003 and .0011. If the snow, 
however, contains as little as 1.6% water content (by weight), the 
loss tangent increases by a factor of 20. Von Hippel (1961) gives 
the loss tangent (at f = 3 GHz) of freshly fallen snow as .0009 
and of hard-packed snow followed by light rain as .0029. 

A typical calculation showing the effect of the snow layer 
on the ice is shown in Figures 15 (a) and 15 (b). If the snow 
is completely dry, it has very little loss and acts as an impedance 
matching layer. For certain frequency and thickness combinations, 
the dependence of the power reflection coefficient on the ice 
thickness completely vanishes. This can be seen for example in 
Figure 15 (a). In this instance, the reflection coefficient shows 
negligible dependence on the ice thickness for frequency of .3GHz 
and a 20 cm snow layer. If the snow, on the other hand, has some 
moisture, it will become more lossy, As the thickness of the 
moist snow increases, the total loss and the emissivity increase^ 
and the ripple disappears. Thus, the variation in the apparent 
brightness temperature as a function of the ice layer thickness 
can be greatly reduced by either dry or wet snow. From Figure 
15 (b)^ it can be seen that for wet snow the apparent brightness 
temperature will also greatly increase with increasing snow layer 
thickness. Even at higher frequencies (33.6 GHz), however, the snow 
cover is often sufficiently transparent that terrain features such 
as field patterns, roads, etc., are all readily identifiable under 
the snow cover, Moore (1974). Thus in practice, unless there is 
a fairly thick cover of wet snow, the variation of brightness with 
ice thickness should remain observable. 

4 . Water on Top of Ice 

During periods of warm weather, some of the surface ice may 
melt forming a layer of water on top of the ice. Changes in 
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Figure 15. Changes in the power reflection coefficient due to 
a snow layer on top of the ice. Note, a quarter 
wavelength in dry snow is .19 meters. 
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power reflection coefficient with increasing water layer thick- 
ness are shown in Figures 16 (a), (b), and (c). It is clear that 
it may be possible to tolerate a millimeter or so of water on top 
of the ice without destroying the Tb versus d curve but not much 
more. The actual power absorbed in the water layer is not large 
as shown in Figures 17 (a), (b), (c), and Figures 18 (a), (b), 
(c), but the power reflection coefficient has changed drastically 
due to changes in the impedance match. The effects of a water 
layer on reflection coefficient at lower frequencies (300 MHz) 
are discussed in Smith (1972). 

5. Slush on Top of Ice 

During weather conditions producing alternate thawing and 
freezing, it is possible that a layer of slush may form on top of 
the ice layer. For lack of published data on the dielectric 
constant for such a layer, recourse must be made to mixture form- 
ulas for slush modeling. 



Figure 16 (a). Changes in the power reflection coefficient due 
to a water layer on top of the ice. 


40 


4-> 

C 

<u 



u • 

•I- cu 

<+- u 

^ *r- 

d) 

O O) 
o jc: 
+-> 
c 

O 4- 

•r- O 
+J 

O Q. 
0) O 

r— 4-> 

4- 

0) C 

t. o 


s- s_ 

O) 0) 

>> 


O fO 


0.r— 


cu &- 
JO. cu 
+-> +-> 
n3 
C 5 


rt3 

tn 

OJ o 
cn+J 
c. 

ca (u 

JO =3 

x_) -a 







Figure 16 (c). Changes in the power reflection coefficient 
due to a water layer on top of the ice. 
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Figure 17 (a-b). Power absorption in the ice/water system as a 
function of aspect angle, assuming no water 
layer on top of the ice. 
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Figure 17 (c). Power absorption in the ice/water system as a 
function of aspect angle, assuming no water 
layer on top of the ice. 


44 



FRACTION OF POWER 


1 . 0 , 



Figures 18 (a-b). Power absorption in the ice/water system as a 

function of aspect angle, assuming a water 
layer on top of the ice. 
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Figure 18 (c). Power absorption in the ice/water system as a 
function of aspect angle, assuming a water 
layer on top of the ice. 


A number of mixture formulas have been derived, but those with 
an adjustable parameter, called the form number, offer the greatest 
flexibility in describing effects such as particle shape, size, 
orientation, distribution, etc.. One such formula is due to Weiner, 
Evans (1965) 


e 

m 


- 1 


+ u 


P(ei - 1) 

+ u “ 


(4-18) 


where 


^m relative complex dielectric constant of the mixture 

e, , g 2 are the relative complex dielectric constants of the 
two mixture components 

p is the fraction of mixture volume occupied by component 1 
u is the form number (o < u < “) 












If the two components are thoroughly mixed, or if one component 
in the form of chunks is suspended in the other (such as small 
pieces of ice floating in waterj^ then u is on the order of 10. 

If the two mixture components are in the form of interleaved layers, 
then the form number takes its limiting value, u approaches zero 
if the_ layer boundaries are oriented orthogonal to the propagation 
direction of the incident field and approaches infinity if the layer 
structure is parallel (i.e., in the latter cas^the field does not 
cross the boundaries between the two constituents, Evans (1965)). 

The relative complex dielectric constants of water and ice 
given in the previous sections were substituted in Equation (4-1 8) 
to calculate possible values of the dielectric constant of slush. 

The real and imaginary parts of the water-ice mixture dielectric 
constant for various form numbers are shown in Figures 19 (a) and 
(b). An equal volume occupancy by the two constituents has been 
assumed in this calculation. The great variation in dielectric 
constant possible for a given volume occupancy and the possibility 
of any fractional composition indicates the difficulty in realistic 
modeling of the actual slush found in the natural state. 

The variation of the power reflection coefficient with increasing 
slush layer thickness was calculated using a representative value 
of p = .5 and u = 10 in the mixture formula. This choice of 
percentage volume occupancy and of form number gives the relative 
complex dielectric constant values listed in Table 2. 

TABLE 2 

Water-Ice Mixture Dielectric Constant Values 


f(GHz) 


Ini(cv^P 

0.1 

13.274 

-0.0406 

0.2 

13.274 

-0.0595 

0.3 

13.274 

-0.0831 

0.4 

13.273 

-0.1079 

0.5 

13.272 

-0.1333 

0.6 

13.272 

-0.1589 


The change in the power reflection coefficient with increasing 
slush layer thickness is shown in Figures 20 (a) and (b). It is 
clear that a few centimeters of slush on top of the ice would 
seriously affect most data processing methods that rely on the 
shape of the power reflection coefficient. It should be noted 
that because the actual properties of the slush layer are not 
known, calculations such as this one can be considered only as 
an indication of the true behavior of the surface. 
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Figures 19 (a-b). The relative complex dielectric constant of 

a water-ice mixture as a function of 
frequency and form number. 
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Figures 20 (a-b). Changes in the power reflection coefficient 

as a function of slush layer thickness on 
top of the ice. 




6. Slush Between Water and Ice 


It may also happen that a mixed layer of water and ice exists 
on the bottom surface of the ice slab, though that is not usually 
the case for a fresh water ice system. The same values as in the 
previous section were used to model the relative complex dielectric 
constant of the slush layer. The power reflection coefficient of 
the ice layer-slush layer-water system were calculated for several 
different slush layer thicknesses and frequencies of 0.4 GHz and 
0.6 GHz. These results are shown in Figures 21 (a) and (b). Two 
different effects can be seen here. At f=0.4 GHz the slush acts 
more like a resonant matching layer which reduces the oscillation 
in the reflection coefficient only at certain thicknesses (4 cm, 
for example). At .6 GHz, there appears to be a steady decline in 
the oscillation of the power reflection coefficient due to 
"smoothing out" of the ice-water boundary. Again, these calcula- 
tions are useful only for indicating qualitative behavior due 
to the large variation possible in slush dielectric constant. 

As stated before, slush layers below the ice are not expected 
in fresh water. 

7. Small Scale Surface Roughness 

Any real ice surface found in the natural state will have 
some surface texture. The effect of this small scale surface 
roughness on a boundary between two media is to cause diffuse 
scattering in all directions at the expense of power in the 
specular direction. The order of power reduction in the 
specular direction is given by. Ruck (1970), 

|R_ |2 ^-{2kh cos6)^ (4-19) 


where 


h is the RMS surface roughness height 

0 is the angle of incidence 

k is the propagation constant in air 

R^. is the reflection coefficient (air- ice) • 

The exponential factor is plotted in Figure 22. Usually the 
criterion h/A .04 defines a "slightly rough surface." Thus, 
for frequencies around 1 GHz for example, the RMS surface rough- 
ness height should be less than 1.2 cm for Equation (4-19) to 
be valid. 
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Figufes 21 (a-b). Changes in the power reflection coefficient as 

a function of slush layer thickness between 
the water and the ice. 
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Figure 22, The 2 reduction of the specular power reflection coefficient 
|R^f , due to slight surface roughness. 

The fraction of incident power scattered diffusely from a slightly 
rough surface may be calculated more precisely using the dif- 
ferential incoherent scattering cross-sections a.. (o»s) and 
a . .( o,s)^ Ruck (1970) and Peake (1971). Consider*^ power density I_ 
iAiident on an area A at an angle 0 .. The total incident power^ 

P. is then ^ 


P. = Iq A cos . (4-20) 

The power density scattered from A at a distance R is 

T ^ ^°.ik ^ .li 

s ^2 

and the total scattered diffuse power is 

2 7T TT^ 2 

Ps = J j R sine^ de^ 

< 1 > 5=0 e ^=0 


(4-21) 


(4-22) 
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and 


P_ = 




2 TT IT/Z 
S S [a 

^3=0 03=0 


jl^ (o,s) + Cjj (o,s)] (4-23) 


sin Qg ®s ^ 


The fraction of power scattered diffusely, is then 
1 


F. = 


d 4it cose Q (jjg-O e^-O jk 


. sin d d 


2 IT ir/2 

In ^ In (o,s) + Qjj (o,s)] 


(4-24) 


For a slightly rough suface, the differential incoherent scattering 
cross-sections are given by Ruck (1970) 


4 . 4 .2 


'hn 


(o,s) = - k h cos 0 . cos 0 a. . 1 I (4-25) 
7 T 1 s nri ' 


°hv 


°vh 


°vv 


4 4 2 2 2 

— k h cos 0 . cos 0 

TT T 


s l^hv 


4 4 2 2 2 , 

7 k h cos 0 . cos 0 ^ 

4 4 2 2 2 

— k h cos 0 . cos 0 ^ 

7 T 1 s ' vv 


(4-26) 

(4-27) 

(4-28) 


where 


- (e - 1) COS <f) 

_ C s 

" ( cos + /77^:^n^“67r(cos 

(4-29) 


- sin (e^ - 1) /e7 - si n^ e 


s ^~c ^ ^c _s 

(cos eV + V^T(e cos' e^ + /e*^"-"Tfn^ e 


s 

(4-30) 
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“vh 


^ . . 2 — 
sin 6 (e^ - 1) /e^ - sin e . 

C C T 

(e COS e. + v^e' ■ sin eTTCcos e + /e™-sin e ) 

C 1 C 1 b b 

(4-31) 


(e(, - 1)(e^ sin sin - cos - s'i e7e^, - sin^ e^) 

'''' (e^ COS e. +/e"^' ^in COS - sin e^) 

(4-32) 

2 2 
I t 

I = ttA ^ e (4-33) 

for a Gaussian correlation coefficient, and 


I = Zvli.^ 


(I + 


(4-34) 


for an Exponential correlation coefficient, where 


t = 




= sin - sin 0^ cos (j)^ 

^ = sin 0^ sin 

s s 

2 7T 

k = is the propagation constant in air 

A 

0. is the angle of incidence 

e^, (j)g are the scattering angles 

h is the RMS surface roughness height 
i is the surface height correlation length 

is the relative complex dielectric constant. 


(4-35) 

(4-36) 

(4-37) 
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PERCENTAGE 


Surfaces having a Gaussian surface height correlation coefficient 
are smoothly curving, with derivatives existing at all points. 
Surfaces with exponential correlation coefficients are jagged 
and have many vertical facets. 

The significance of these fractional power calculations for 
the perceived brightness of the surface may be seen by considering 
the total reflection from the surface in two parts. The fraction 
F. of the incident power that is scattered diffusely above the 
ice surface is given in Equation (4-24). Numerical integration 
of this expression for both polarizations gives the results shown in 
Figures 23 (a), (b), {c\ and (d). It should be emphasized that 
these curves indicate only the amount of power diffusely scattered 
into the half-space above the ice surface. 



ANGLE OF INCIDENCE (DEGREES ) 


Figure 23 (a). The percentage of the incident power that is 
scattered diffusely above the slightly rough 
surface. 
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PERCENTAGE 




Figures 23 (b-c). The percentage of the incident power that is 

scattered diffusely above the slightly rough 
surface. 
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DIFFUSE SCATTERING -ICE 


h = 0.04X 



Figure 23 (d). The percentage of the incident power that is 
scattered diffusely above the slightly rough 
surface. 


The actual change in apparent brightness temperature caused 
by the surface being slightly rough may then be computed using 
results of Valenzuela (1970). He showed that the effective re- 
flection coefficients for horizontal and vertical polarization are 




k^.t( Ec ‘ 




— 2^2 
- T ) 



(4-38) 


(4-39) 


where 
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(4-40) 




k (s^ - l)(e^ 


- 



where 

is the relative complex dielectric constant 
w(p,q) is the spectral density of the surface profile 
k is the free space propagation constant, 
and 

^ = c os e^. T = sin e. 

b = /k ^ C =yk ^ 

d = b + c D = be + p^ + 

r H 
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The change in the apparent brightness temperature due to the 
surface roughness is then 


= (iRl^ - |R'I^ - Fj) T. 


(4-44) 


where 


^Tb = T|^ (slightly rough surface) - Tj^ (flat surface). 


Because the sky brightness is quite small for the frequencies of 
interest here, the change in brightness due to reflected components 
has been assumed to be negligible in Equation (4-44). Numerical 
integrations of Equations (4-42) and (4-43) were performed to 
calculate at. . Representative results are shown in Figures 24 (a), 
(b), (c), and (d). These results can be considered only indicative 
of the brightness changes to be expected due to slight roughness 
on the surface. As can be seen from the figures, the slightly rough . 
surface will look a few degrees warmer, than the smooth surfaces 
for scan angles of 30®, or less. 


The rather complicated reflection coefficients R/, R' for a 
slightly rough surface are often approximated by the BimpYer expres- 
sion (4-1 9j where the exponential term is the power loss factor caused 
by roughness for a perfectly conducting surface. That the approxi- 
mation is quite a good one may be judged from Table 3, which displays 
the correct and approximate reflection coefficients for an ice 
surface with h/X = .04 for several values of correlation length 


TABLE 3 

Power Reflection Coefficient of Slightly Rough Surface 
Gaussian Correlation Coefficient - Horizontal Polarization 


-(2kh cose)^ 

9 IRJ^ iRj'e 


0° 

0.0804 

0.0625 

10° 

0.0832 

0.0651 

20° 

0.0921 

0.0737 

30° 

0.1091 

0.0903 

40° 

0.1382 

0.1191 

50° 

0.1864 

0.1679 

60° 

0.2668 

0.2505 

70° 

0.4016 

0.3899 

80° 

0.6274 

0.6626 

90° 

1.0000 

1.0000 


&k=1.0 


&k=4.0 


K \^ 

|r;,i^ 

0.0588 

0.0626 

0.0614 

0.0613 

0.0654 

0.0642 

0.0693 

0.0744 

0.0732 

0.0847 

0.0917 

0.0907 

0.1115 

0.1211 

0.1210 

0.1572 

0.1699 

0.1717 

0.2352 

0.2511 

0.2557 

0.3699 

0.3876 

0.3949 

0.6025 

0.6173 

0.6246 

1.0000 

1.0000 

1.0000 


59 


Gaussian Correlation Coefficient - Vertical Polarization 




-(2kh cos0)^ 

Ak=1.0 

Ak=2.0 

aM.O 

6 

|R 1^ 

iRvl^e 

|R' 1^ 
1 v‘ 



0° 

0.0804 

0.0625 

0.0588 

0.0626 

0.0614 

10° 

0.0777 

0.0608 

0.0572 

0.0609 

0.0598 

20° 

0.0694 

0.0555 

0.0522 

0.0557 

0.0547 

30° 

0.0554 

0.0459 

0.0433 

0.0462 

0.0453 

40° 

0.0363 

0.0313 

0.0299 

0.0318 

0.0309 

50° 

0.0145 

0.0131 

0.0135 

0.0137 

0.0129 

60° 

0.0001 

0.0001 

0.0023 

0.0004 

0.0001 

70° 

0.0266 

0.0258 

0.0275 

0.0232 

0.0252 

80° 

0.2091 

0.2075 

0.2034 

0.1981 

0.2040 

90° 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 


Exponential 

Correlation Coefficient - Horizontal Polarization 



-(2kh cose)^ 

')ik=1.0 

Ak=2.0 

Jtk=4.0 

0 

iRhl^ 

iRhl^e 


IR^I^ 

IR^I^ 

0° 

0.0804 

0.0625 

0.0570 

0.0597 

0.0609 

10° 

0.0832 

0.0651 

0.0595 

0.0623 

0.0636 

20° 

0.0921 

0.0737 

0.0676 

0.0709 

0.0724 

30° 

0.01091 

0.0903 

0.0832 

0.0874 

0.0892 

40° 

0.1382 

0.1191 

0.1103 

0.1159 

0.1183 

50° 

0.1864 

0.1679 

0.1564 

0.1638 

0.1675 

60° 

0.2668 

0.2505 

0.2348 

0.2446 

0.2497 

70° 

0.4016 

0.3899 

0.3699 

0.3811 

0.3876 

80° 

0.6274 

0.6226 

0.6028 

0.6124 

0.6182 

90° 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 


Exponential 

Correlation Coefficient - Vertical Polarization 



“(2kh cose)^ 

JII<=1.0 

Jlk=2.0 

aM.o 


In |2 

Ir^ |2 

1 _ I 1 2 

IrM |2 

Ir.. l2 

e 


R e 

R' 

R' 

R' 


1 V ‘ 

1 V ' 

' v‘ 

1 V ' 

1 v' 

0° 

0.0804 

0.0625 

0.0570 

0.0597 

0.0609 

10° 

0.0777 

0.0608 

0.0555 

0.0581 

0.0593 

20° 

0.0694 

0.0555 

0.0508 

0.0530 

0.0542 

30° 

0.0554 

0.0459 

0.0426 

0.0439 

0.0448 

40° 

0.0363 

0.0313 

0.0309 

0.0302 

0.0305 

50° 

0.0145 

0.0131 

0.0176 

0.0134 

0.0128 

60° 

0.0001 

0.0001 

0.0108 

0.0020 

0.0004 

70° 

0.0266 

0.0258 

0.0405 

0.0279 

0.0258 

80° 

0.2091 

0.2075 

0.2161 

0.2054 

0.2048 

90° 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 
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Figures 24 (a-b). Changes in the apparent brightness temperature 

due to slight roughness on the surface. Reflected 
components have been neglected. 
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24 (c-d). Changes in apparent brightness temperature due 
to slight roughness on the surface. Reflected 
components have been neglected. 


8. Multilayerinq of Ice - 

Uarge Scale Roughness^ 

The small scale roughness theory of the preceding section can 
predict the effect of surface roughness between two media on the 
order of a few centimeters or less (for the frequencies concerned 
in this report). The actual lake ice could very well have surface 
variations much greater than this if, for example, any rafting is 
present. The radiometer antenna could also see a continuum of 
ice layer thicknesses if it was observing the ice close to open 

water. The effect of the presence of a number of different 
thicknesses within the antenna footprint on the apparent bright- 
ness temperature is shown in Figure 25. The "9.5 cm" curve, 
for example, shows the temperature the antenna would see if on 
top of any given ice layer there were additional slabs ranging 
in thickness from 0.5 to 9.5 cm, with equal weighting as seen by 
the antenna. The presence of multiple thicknesses tends to smooth 
out the oscillation in brightness. The increase in the effective 
total thickness is reflected in the lateral shift of the sinusoidal 
curve. The effect of the smoothing in the frequency domain is shown 
in Figures 26 (a) and (b). The separation of the nulls along the 
frequency axis is inversely proportional to the spread in thickness 
as predicted by Equation (2-4). For example, in Figure 26 (b) when 
the thickness spread increases from 4 cm to 10 cm, the location of 
the first null changes from 2.0 GHz to 0.8 GHz. Furthermore, the 
null locations depend only on the multi layering, and not on the basic 
ice thickness, so that a 10 cm thickness variation places a null 
around 0.8 GHz regardless of the average ice thickness. It is clear 
that the presence of multiple thicknesses will have an adverse effect 
on most ice thickness prediction methods that depend on the presence 
of the interference phenomenon, and place an upper limit on the 
maximum useful frequency, f - c/(3/^ Ad) where Ad is the range 
of thicknesses within the beam. 

9. Air Bubbles in Ice 


The ice model used in the previous calculation is quite low 
in loss, and applies to ice free of cracks, bubbles, inclusions, 
etc.. Unfortunately, that is not always the case for ice found in 
the natural state. Ice samples taken from the Great Lakes, for 
example, can show a high air bubble content. Measurement, Vickers 
(1975), of the relative complex dielectric constant of these 
samples produces real parts that vary between 3.2 and 3.0 with 
losses from 1 dB/m to approximately 4 dB/m at 1 GHz. This loss 
variation corresponds to the imaginary part of the dielectric 
constant varying from -.02/f to -.08/f (f is in GHz). The classic 
mixture formulae correctly predict the change in the real part of 
the dielectric constant but not in the imaginary part, and 
therefore cannot be used in modeling air bubbles in ice. The 
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Changes in the apparent brightness temperature d.r 
to presence of additional ice slabs on top of the 
basic ice layer. 
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Figure 26 (a). Changes in the apparent brightness tempers. ture 
due to presence of ice multi-thickness. 
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Figure 26 (b). Changes in the apparent brightness temperature due 
to presence of ice multi- thickness. 


effect of the increased lossiness on the apparent brightness temper- 
ature for two different ice thicknesses as a function of frequency 
is shown in Figures 27 (a) and (b). Typical effects of lossy ice 
and multi layering on the apparent brightness temperature at .4 GHz 
is shown in Figure 28. The "harmonic content" changes and the 
'‘dc" level shifts considerably, but the period of oscillations re- 
mains about the same. Any data processing method will have to ad- 
just to this fact if it is to cope with lossiness and multi layering. 


65 


800 



Figure 27 (a). Changes in the apparent brightness temperature due to 
presence of air inclusions in ice. 
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Figure 27 (b) . 


Changes in the apparent brightness temperature due 
to presence of air inclusions in ice. 
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Figure 28, The apparent brightness variation with increasing ice loss 
and ice multi layering. 


CHAPTER V 


DATA PROCESSING FOR CLEAR SINGLE LAYER ICE 
0. Introduction 


In this chaptei} we consider two types of realizations of 
automatic ice thickness measuring schemes. The first is based 
on measurements at frequencies low enough (say 150 MHz or less) 
so that the apparent brightness temperature is dominated by the 
reflected cosmic radiation. At such low frequencies, for ice 
up to a meter or so thick, only two or three frequency channels 
are needed for an unambiguous thickness determination. The 
second type of system, operating at frequencies where the 
brightness temperature is dominated by emission (about 300 MHz 
or higher), requires in its most direct form, N = d /d . +1 

channels. The object of the work described here, hoWi^er^^*^ 
is to develop data processing schemes which will permit unam- 
biguous thickness determination using the least possible number 
of frequency channels consistent with a reasonable constraint 
on the probable ice thickness error. All calculations in this 
chapter will assume the ice to be clear and of uniform thickness. 

1 . Low Frequency Measurements 

This section is basically concerned with how to measure ice 
as much as a meter thick, using as few frequencies as possible. 

It will be seen that the requirement for minimum number of 
frequency channels forces the frequencies to be very low. At low 
frequencies, however , the galactic noise temperature is a very 
large, imprecisely known quantity, and the data processing 
methods have to be designed to cope with this. Although in 
principle, by working on the first half cycle of the reflection 
coefficient vs. ice thickness curve, only a single frequency would 
suffice, we do not consider such systems here because at the 
frequencies required (less than 50 MHz), problems of antenna 
design and cosmic noise uncertainty make them impractical. 

A. Two Antenna System 

From Figure 6, it can be seen that unambiguous resolution of 
ice thicknesses up to one meter, by means of brightness measure- 
ments at two or three frequencies^ requires that the frequencies 
be on the order of 0.1 GHz. In this frequency range^ the galactic 
noise temperature has great variation but is in general very 
large as shown in Figures 4 and 5. The actual value of the 
galactic noise temperature depends on the position in the sky. 
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the antenna beamwidth, etc. It is on the order of many thousands 
of degrees near the galactic center, and it falls o'f'f to many 
hundreds of degrees in the direction of the galactic pole. The 
physical temperature of ice, on the other hand, is 273° K or less, 
and so the reflection mechanism clearly dominates the emission 
at these low frequencies. The apparent brightness temperature 
was given as 


= (1 - |R|^) T. + |Rl^ (Tg + Tj,) , (3-1) 

SO that if T. 

g a 1 

then iRl^ (T^ + T^) . (5-1) 

It should be noted that the atmospheric noise temperature 
contribution T is on the order of few degrees K, and thus 
may be considered at these low frequencies negligible in 
comparison with T . The apparent brightness temperature can 
be viewed then as^the power reflection coefficient multiplied 
by the large but variable galactic noise. Hence, it would 
be desirable if the data processing method did not have to 
rely on accurate, absolute temperature measurements. 

One obvious approach to the problem is to measure both 
T + T , and T. . Thus, T. is measured with a downward 
l8oking antenna, and the^temperature incident upon the ice 
layer from above, namely T » is measured with an upward 
looking antenna. As can be seen from Equation (5-1), the ratio 
of T. to T + T will be an approximation to the power re- 
flection cSeffiSient, i.e. 

1R|^ = {Tj, - ip/(Tg + - Tp^ydg + y . (5.2) 

A slightly better approximation would be to assume that the 
unknown ice temperature T^ is 270° K, so that 

iRl^^ (T,^ - 270)/(Tg + Tg - 270) . (5-3) 

A theoretical, or pre-calculated, power reflection coefficient 
curve can then be used as a training set for identification 
of each measured unknown ice thickness. Thus, the training 
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2 

set for a single frequency realization is just a graph of |R | 
vs thickness. For example, (see Figure 13) at 0.1 GHz^only a 
small range of thicknesses (0cm -40cm) can be measured unambig- 
uously. To extend the unambiguous range, two frequencies may 
be used; the training set can then be represented as a plane 
curve with ice layer thickness as a parameter. For example. 
Figure 29 (a) is the locus of 199 two-dimensional vector tips 
representing thicknesses from 1 to 100 cm, in 0.5 cm steps, 
for two frequencies of 0.9 GHz and 0.11 GHz. By using two 
frequencies, the unambiguous thickness measurement range is 
now approximately 15 cm to 70 cm. To implement such a system, 
one must then, at each of the two frequencies, measure the 
brightness seen by an upward looking antenna and a downward 
looking antenna, and use the up/down brightness ratios as 
estimates of the reflection coefficients, these being the two 
components of the measurement vector. 



Figure 29 (a). Pre-computed power reflection coefficients at 
two frequencies to be used as a training set 
for ice thickness prediction. 
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The actual thickness estimate is accomplished using a 
minimum distance algorithm, Tou and Gonzales (1974), (often 
known as the nearest neighbor rule). The measurement vector 
from an unknown ice thickness is compared to each of the 
vectors in the training set to determine the best fit. The 
criterion for best fit is simply the minimum geometric distance 
between the measurement vector and each of the other vectors 
in the training set. Once the minimum distance has been 
found, the identification process is complete, and the measured 
ice thickness is assumed to be that of the vector in the 
training set which produced the best fit. 

As mentioned above, even for frequencies as low as 0,09 
GHz and 0.11 GHz an ambiguous crossover point exists in the 
training set. For example, (see Figure 29 (a)) a measurement 
vector from ice 10 cm thick could easily be interpreted as 
ice 84 cm thick. Clearly, the set of ambiguous thicknesses 
will depend on the measurement error. To reduce the ambiguity 
range, at the expense of system complexity, one must include 
data at a third frequency, thus, increasing the dimension of the 
measurement space from two to three. The locus, or training 
set_,must now be plotted in a three dimensional space. The 
two dimensional locus (e.g. Figure 29 (a)) would then be the 
projection of this three dimensional curve onto a plane. Be- 
cause the three dimensional locus is difficult to present (two 
projections of such a locus are shown in Figures 29 (b) and (c)) 
we give here only the results for some representative system 
realizations. A typical performance result for temperature 
measurements 2at three frequencies and for a galactic temperature 
of T - 20 /t is shown in Figure 30. Depicted is the predicted 
vs actual ice thickness using the above outlined method with 
a 199 member training set and with the measured power reflection 
coefficient based on the approximation [R f = T./(T + T ). The 
average error of 1.475 cm is the sum of deviatTons^from^the 45 
degree line divided by the number of thicknesses. Clearly, 
apart from a small quantization error, the minimum distance 
algorithm should give almost perfect results if the exact 
expression. Equation (5-2), were used. However, the error curve 
shown in Figure 30 can be considered typical of the results that 
might be achieved with such a system for systematic errors of 
order of 300K° in measuring the upward and downward brightness 
temperature. 

As pointed out before, this data processing method depends 
on having T >> T.. Thus, as the galactic noise temperature be- 
comes smaller, thJ performance will deteriorate. The effect is 
shown in Figure 31, which repeats the "average error" calculation 
of Figure 30 for a range of galactic temperatures. When G has 
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Figure 30. Ice thickness prediction errors using low frequency 

two-antenna measurement system for galactic temperature 
factor G=20. 
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AVERAGE ERROR (cm) 



GALACTIC TEMPERATURE FACTOR G 

rfgure 37. Ice thickness prediction errors using low frequency two-antenna 
measurement system as a function galactic temperature factor G. 




decreased to about 10, this data processing method fails at 
frequencies near 0.1 GHz for the kind of measurement errors 
represented by the approximation in Equation (5-2). In practice, 
one would use the better approximation, i.e. assume T. - 270° K, 
but would be subject to both random and systematic measurement 
errors. Furthermore, since the galactic temperature is a 
function of position in the sky, if the downward looking antenna 
is to scan, then so must the upward looking antenna. Clearly, 
a far more practical arrangement would be one which requires 
only a single downward looking antenna. 

B. One Antenna System 

If the airborne measurement system has only one downward 
looking antenna, then it is not possible to measure the galactic 
plus atmospheric contributions to the apparent brightness temper- 
ature explicitly. The effect of the unknown and variable 
galactic noise temperature can be minimized by a data processing 
approach that utilizes the brightness temperature ratio at a 
pair of frequencies as the component of the measurement vector. 
Thus, from three frequency channels, one may find two temperature 
ratios to obtain a two component measurement vector. Assuming 
again, that the frequency is low enough so that T is very large 
compared to T. and T can be neglected, then the apparent 
brightness terliperatures measured at three frequencies can be 


considered to be 


^ [R {fpd)l^ G/f^ 

(5-4) 


(5-5) 

T^(f3) IR (fa.d) 1^ G/f3 

(5-6) 

2 

where [R (f ,d) | is the power reflection coefficient at frequency 
f and the ice layer thickness d. Expressions for the power 
reflection coefficient at either polarization are given in 
Equation ( 4 - 7 ) and Equation (4-12). The two components of the 
measurement vector then become for example. 

vf^) jR(fi.d)|2 

|R(f3,d)|^ VT/ 

(5-7) 


and 
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2.7 


(5-8) 



|2 

N 

lR(f3.d)i 

7 ' 

i^2/ 


To the extent that T. >> T. » T , the two numbers represented 
by these two equations areVunctfons of ice layer thickness but 
not of the absolute value of the galactic noise temperature. 

The invariance of the measured temperature ratios with 
changes in the galactic temperature can also be seen from a 
geometric consideration. If a plot is made of T. (f=0.09 GHz,d) 
vs T. (f=0.11 GHZjd) for example, for all ice thicknesses d, the 
curve will be similar in form to the curve of Figure 29 (a), 
since at low frequencies the apparent brightness temperature is 
approximately proportional to the power reflection coefficient. 

If the galactic temperature factor G is increased, both temper- 
atures at each ice thickness would increase proportionally, 
and the shape of the curve would remain the same. Hence, the 
slope of any vector from the origin to any given ice thickness on 
the temperature curve would remain the same. The value of the 
slope is just given by Equation (5-7). This approximate in- 
variance of the measured temperature noise ratios is the basts 
for the data processing method of this section. 

Each of the vectors in the training set is constructed 
using a theoretical temperature ratio for each vector component. 

A general expression for the two-dimensional training vector as 
a function of ice layer thickness d is 



where 

S 0 - |R(f^.d)l^)T. + |R(fi,d)|^(T^(G.f^) 

^ 0 - |R(f3,d)|^)T. + |R(f3,d)l^(Tg(G,f3) + Tg(f3)) 

(5-9) 


and 
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In the data processing method of the preceding section, the 
training set was formed using theoretical estimates of the power 
reflection coefficient. In this case, an assumed value of G 
must be incorporated in the training set. One such training 
set is shown in Figure 32. This, as before, is the locus of 
199 vectors representing ice layer thickness from 1 cm to 100 cm 
in 0.5 cm increments. One component of each vector is given 



Figure 32. Ratios of pre-computed apparent brightness temperatures 

to be used as a training set for ice thickness prediction. 


by S-j , the other component by Sp. Here the maximum possible 
galactic noise (G=40) has been used for calculating the vectors 
of the training set. The measurements are simulated by computing 
an "observed" brightness temperature using a G different from 
the one used in calculating the training set. The measurement 
vectors are then constructed in the same way as the training set 
vectors. The remainder of the data processing uses the minimum- 
distance algorithm as outlined before. Once a two-dimensional 
measurement vector is formed from three temperature measurements, 
the measurement vector is compared to each of the 199 vectors 
in the training set to find the best fit, the criterion again 
being the minimum geometric distance between the measurement 
vector and each of the training set vectors. 

If the apparent brightness temperature measurements are 
made in an environment having the same galactic temperature as 
assumed in the construction of the training set, the performance 
of this data processing method should be errorless, in principle 
If, on the other hand, the actual galactic noise temperature 
differs from that assumed in the training set, the performance 
of the data processing method will deteriorate. For example, 
if the training set of Figure 32, which was constructed assuming 
the maximum galactic noise (G=40), is used to identify measure- 
ment vectors that have been made in an environment having an 
actual galactic noise temperature T =20 /t i.e., G=20, the 

performance results are shown in Figure 33. The average error 
of 0.832 cm is^ as before, the sum of the deviations of the 
predicted thickness curve from the 45 degree line divided by 
the total number of thicknesses (199). 

It is clear that the closer the galactic noise factor G 
used in construction of the training set is to the G assumed 
to be in the measurement environment, the better will be the 
performance of this data processing method. This variation in 
performance is shown in Figure 34. Plotted here are the average 
errors in predicted thickness for training sets (constructed with 
G = 40, = 30, =20) as a function of the actual galactic temper- 
ature. 

Comparing Figure 31 and Figure 34, it can be seen that both 
methods are comparable in performance, but the latter method has 
the advantage of requiring only one antenna. In both methods, 
there is no need for accurate absolute temperature measurements 
since calibration errors have only a small effect on the bright- 
ness ratios used to form the measurement vector. Both methods, 
however, do depend upon the galactic noise temperature being much 
larger than the physical temperature of the ice, and the performance 
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ACTUAL ICE THICKNESS ( Cm ) 

Figure 33. Ice thickness prediction errors using low frequency 

one-antenna measurement system; training set calculated 
using G = 40, "measurements" made in an environment 
having G = 20. 
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Figure 34. Ice thickness prediction errors using low frequency 

one-antenna measurement system; training sets calculated 
using G = 40, G = 30, G = 20, "measurements" made in 
environment having all possible values of G {2^G:i40). 
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of either method deteriorates rapidly if the galactic temperature 
factor G decreases below 10 for frequencies around 0,1 GHz. 

The design of a suitable airborne antenna for frequencies this 
low is a formidable problem which can be greatly alleviated by 
going to higher frequencies. At such frequencies (Say f > 300 MHz), 
the galactic noise temperature becomes very small, and the apparent 
brightness temperature is dominated by the emission from the water - 
ice layer system. 

2. High Frequency Measurements 

It was shown in the preceding section that ice up to a meter 
thick could be measured in principle with two or three frequencies 
in the 100 MHz region. At these frequencies, however, one must 
accommodate a difficult antenna design and a variable cosmic noise 
level. Both of these drawbacks are circumvented by going to higher 
frequencies (up to 1.0 GHz). This section will be concerned with 
data processing methods for temperature measurements made at these 
higher frequencies. 

A. Minimum - Distance Pattern Classification 

From Figure 6 it is seen that for frequencies on the order 
of 0.1 GHz the power reflection coefficient, and hence the apparent 
brightness temperature, experiences approximately one periodic 
oscillation over the span of ice layer thicknesses ranging from 1 
cm to 100 cm. If the power reflection coefficient of the frequencies 
on the order of 1.0 GHz was plotted in the manner of Figure 29 (a), 
it is clear that one would obtain a Lissajous type figure with many 
crossover points, giving rise to many ambiguities. Use of such a 
training set for identifying two or three dimensional measurement 
vectors would not be feasible. It was noted earlier, however, that 
the addition of more dimensions to the measurement space tends to 
"untie" the intersections in the training set. To be investigated 
here,are the number of frequencies required for a satisfactory 
thickness estimate consistent with a realistic measurement error, 
the size of the frequency sample span, the frequency sample dis- 
tribution within the span, and the location of the sample span 
within the frequency spectrum. 

The training sets are calculated using the theoretical model 
for the apparent brightness temperature, i.e. Equation (3-1). 

The minimum value of the galactic temperature (G=2) has been 
assumed in the training set calculations. (At frequencies of 1,0 
GHz or higher^ the galactic noise makes an insignificant con- 
tribution to the brightness temperature of the ice). 
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In order to simulate systematic errors in the measuring 
system due to ignorance of such parameters as the actual ice 
temperature, the exact complex dielectric constant of ice, the 
effects of small scale surface roughness, uncorrected instrument 
drift etc., a constant temperature error, was added to the 
theoretical temperatures at every frequency. The measurement 
vectors were then constructed from these "corrupted" theoretical 
temperature values. The multi -dimensional training set vector is 
then 


(5-11) 


where the brightness temperature S, for frequencies f^ is 

= (I - lR(fk.d)l^)T.+ |R(fk.d)|^(Tg(G.f^) + y^) . 

(5-12) 

An arbitrary measurement vector is then represented by 


(5-13) 


where 

+ AT . (5-14) 

The minimum distance algorithm is applied just as in the data 
processing methods of the low frequency measurement section. 
Each of 199 measurement vectors representing the ice thickness 
of interest is compared to each of the 199 vectors in the 
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TABLE 4 


Minimum-Distance Pattern Classifiction 
(5 K“ Constant Temperature Error) 

Frequency (GHz) Average ice 

Thickness 
Error (cm) 


1.00 

1.09 










29.198 

1.00 

1.05 

1.09 









5.269 

1.00 

1.03 

1.06 

1.09 








2.003 

1.00 

1.02 

1.04 

1.06 

1.09 







1.962 

1.00 

1.01 

1.02 

1.04 

1.06 

1.09 






2.045 

1.00 

1.01 

1.02 

1.03 

1.04 

1.06 

1 

.09 




2.289 

1.00 

1.01 

1.02 

1.03 

1.04 

1.05 

1 

.06 

1.09 



2.289 

1.00 

1.01 

1.02 

1.03 

1.04 

1.05 

1 

.06 

1.07 

1.09 


2.442 

1.00 

1.01 

1.02 

1.03 

1.04 

1.05 

1 

.06 

1.07 

1.08 

1.09 

2.402 

1.00 

1.18 










30.073 

1.00 

1.10 

1.18 









9.178 

1.00 

1.06 

1.12 

1.18 








2.875 

1.00 

1.04 

1.08 

1.12 

1.18 







0.148 

1.00 

1.02 

1.04 

1.08 

1.12 

1.18 






0.221 

1.00 

1.02 

1.04 

1.06 

1.08 

1.12 

1 

.18 




0.181 

1.00 

1.02 

1.04 

1.06 

1.08 

1.10 

1 

.12 

1.18 



0.148 

1.00 

1.02 

1.04 

1.06 

1.08 

1.10 

1 

.12 

1.14 

1.18 


0.083 

1.00 

1.02 

1.04 

1.06 

1.08 

1.10 

1 

.12 

1.14 

1.16 

1.18 

0.113 

1.00 

1.36 










30.957 

1.00 

1.20 

1.36 









10.372 

1.00 

1.12 

1.24 

1.36 








3.178 

1.00 

1.08 

1.16 

1.24 

1.36 







0.246 

1.00 

1.04 

1.08 

1.16 

1.24 

1.36 






0.000 

1.00 

1.04 

1.08 

1.12 

1.16 

1.24 

1 

.36 




0.000 

1.00 

1.04 

1.08 

1.12 

1.16 

1.20 

1 

.24 

1.36 



0.000 

1.00 

1.04 

1.08 

1.12 

1.16 

1.20 

1 

.24 

1.28 

1.36 


0.000 

1.00 

1.04 

1.08 

1.12 

1.16 

1.20 

1 

.24 

1.28 

1.32 

1.36 

0.000 

1.00 

1.63 










29.186 

1.00 

1.35 

1.63 









9.221 

1.00 

1.21 

1.42 

1.63 








3.877 

1.00 

1.14 

1.28 

1.42 

1.63 







0.786 

1.00 

1.07 

1.14 

1.28 

1.42 

1.63 






0.000 

1.00 

1.07 

1.14 

1.21 

1.28 

1.42 

1 

.63 




0.000 

1.00 

1.07 

1.14 

1.21 

1.28 

1.35 

1 

.42 

1.63 



0.000 

1.00 

1.07 

1.14 

1.21 

1.28 

1.35 

1 

.42 

1.49 

1.63 


0.000 

1.00 

1.07 

1.14 

1.21 

1.28 

1.35 

1 

.42 

1.49 

1.56 

1.63 

0.000 
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TABLE 4 (continued) 


Frequency (GHz) Average ice 

Thickness " 
Error (cm) 


1, 

.00 

1.90 








23.950 

1 

.00 

1.50 

1.90 







12.090 

1 

.00 

1.30 

1.60 

1.90 






4.955 

1 

.00 

1.20 

1.40 

1.60 

1.90 





3.628 

1 

.00 

1.10 

1.20 

1.40 

1.60 

1.90 




2.580 

1 

.00 

1.10 

1.20 

1.30 

1.40 

1.60 

1.90 



1.987 

1 

.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.90 


1.834 

1 

.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.90 

1.834 

1 

.00 

1.10 

1.20 

1.30 

1.40 

1.50 

1.60 

1.70 

1.80 

1.90 1.849 
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TABLE 5 


Mini mum- Distance Pattern Classification 
(5 K° Constant Temperature Error) 


Frequency (GHz) 


3.00 

3.09 





3.00 

3.05 

3.09 




3.00 

3.03 

3.06 

3.09 



3.00 

3.02 

3.04 

3.06 

3.09 


3.00 

3.01 

3.02 

3.04 

3.06 

3.09 

3.00 

3.01 

3.02 

3.03 

3.04 

3.06 

3.00 

3.01 

3.02 

3.03 

3.04 

3.05 

3.00 

3.01 

3.02 

3.03 

3.04 

3.05 

3.00 

3.01 

3.02 

3.02 

3.04 

3.05 

3.00 

3.18 





3.00 

3.10 

3.18 




3.00 

3.06 

3.12 

3.18 



3.00 

3.04 

3.08 

3.12 

3.18 


3.00 

3.02 

3.04 

3.08 

3.12 

3.18 

3.00 

3.02 

3.04 

3.06 

3.08 

3.12 

3.00 

3.02 

3.04 

3.06 

3.08 

3.10 

3.00 

3.02 

3.04 

3.06 

3.08 

3.10 

3.00 

3.02 

3.04 

3.06 

3.08 

3.10 

3.00 

3.36 





3.00 

3.20 

3.36 




3.00 

3.12 

3.24 

3.36 



3.00 

3.08 

3.16 

3.24 

3.36 


3.00 

3.04 

3.08 

3.16 

3.24 

3.36 

3.00 

3.04 

3.08 

3.12 

3.16 

3.24 

3.00 

3.04 

3.08 

3.12 

3.16 

3.20 

3.00 

3.04 

3.08 

3.12 

3.16 

3.20 

3.00 

3.04 

3.08 

3.12 

3.16 

3.20 

3.00 

3.63 





3.00 

3.35 

3.63 




3.00 

3.21 

3.42 

3.63 



3.00 

3.14 

3.28 

3.42 

3.63 


3.00 

3.07 

3.14 

3.28 

3.42 

3.63 

3.00 

3.07 

3.14 

3.21 

3.28 

3.42 

3.00 

3.07 

3.14 

3.21 

3.28 

3.35 

3.00 

3.07 

3.14 

3.21 

3.28 

3.35 

3.00 

3.07 

3.14 

3.21 

3.28 

3.35 


Average ice 
Thickness 
Error (cm) 

28.151 

4.995 

2.899 

3.148 

3.425 


3.09 3.452 

3.06 3.09 3.480 

3.06 3.07 3.09 3.342 

3.06 3.07 3.08 3.09 3.286 


28.749 

10.754 

8.706 

0.138 

0.221 

3.18 0.206 

3.12 3.18 0.178 

3.12 3.14 3.18 0.095 

3.12 3.14 3.16 3.18 0.083 

24.771 

10.209 

10.530 

0.050 


0.013 

3.36 0.013 

3.24 3.36 0.013 

3.24 3.28 3.36 0.013 

3.24 3.28 3.32 3.36 0.050 


33.844 

15.013 


4.352 

0.372 

0.307 

3.63 0.402 

3.42 3.63 0.000 

3.42 3.49 3.63 0.025 

3.42 3.49 3.56 3.63 0.000 
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TABLE 5 (continued) 


Average ice 
Thickness 
Error (cm) 


Frequency (GHz) 


3.00 

3.90 



3.00 

3.50 

3.90 


3.00 

3.30 

3.60 

3.90 

3.00 

3.20 

3.40 

3.60 

3.00 

3.10 

3.20 

3.40 

3.00 

3.10 

3.20 

3,30 

3.00 

3.10 

3.20 

3.30 

3.00 

3.10 

3.20 

3.30 

3.00 

3.10 

3.20 

3.30 


3.90 

3.60 

3.40 

3.90 

3.60 

3.90 


3.40 

3.50 

3.60 

3.90 

3.40 

3.50 

3.60 

3.70 

3.40 

3.50 

3.60 

3.70 


28.407 

10.822 

11.538 

0.997 

0.435 

0.435 

0.435 

3.90 0.435 

3.80 3.90 0.435 
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training set to yield a performance graph similar in form to 
Figures 28, 30, 33. In this case, only the average ice thickness 
prediction errors are tabulated here. The results are shown in 
Table 4 and Table 5. As can be seen from the tables, temperature 
measurements for an assumed constant temperature error of 5K° 
at five or more frequencies spanning approximately 0.5 GHz provide 
good results. A slight deterioration in performance is seen as 
the sample span is moved to higher frequencies. This is to be 
expected since the vector locus in the training set has an 
increasing number of twists, turns, and overlaps. 

The cal culationsj whose results are shown in Table 4 and 
Table 5, were made assuming a constant measurement error at of 
5 K°. As this error increases, it can be expected that the average 
ice thickness prediction error will also increase. This rate of 
increase for two different sets of measurements is shown in Table 6. 
The average ice thickness prediction errors are given for increasing 
mi form positive errors, increasing uniform negative errors, and 
increasing errors that alternate in sign from one measurement 
vector component to another for only the nine frequency measurement 
set. From these results^ it can be seen that the performance in 
the face of increasing measurement errors is much better for the 
nine frequency system than for the system that utilizes five 
temperature measurements in the identification process. 

TABLE 6 

Minimum-Distance Pattern Classification 

Variable Temperature Measurement Errors For A Five Frequency 
And Nine Frequency System (Normal Incidence) 


Measurement 

Average Ice Thickness 

(Error in cm) 

Error AT(°K) 

2.00(.06)2.24 GHz 

2.00(. 

03)2.24 GHz 



AT 

AT(-l)^* 

0.0 

0.000 

0.000 

0.000 

2.0 

0.000 

0.000 

0.000 

4.0 

0.040 

0.020 

0.000 

6.0 

0.417 

0.299 

0.000 

8.0 

1.497 

0.580 

0.000 

10.0 

2.643 

0.882 

0.020 

15.0 

5.399 

1.709 

0.630 

20.0 

9.020 

2.633 

1.040 


* Constant temperature error alternating in sign for each component 
frequency. 
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TABLE 6 (continued) 

Error AT(°K) 2.00( .06)2. 24 GHz 2. 00( .03)2.24 GHz 




AT 

AT(-l) 

-2.0 

0.000 

0.000 

0.000 

-4.0 

0.040 

0.020 

0.000 

-6.0 

0.437 

0.281 

0.000 

-8.0 

1.560 

0.603 

0.000 

-10.0 

2.872 

0.899 

0.020 

-15.0 

6.158 

1.641 

0.640 

-20.0 

8.962 

2.332 

1.060 


Clearl)^ when the errors in AT are systematic (of uniform sign for 
all of the measurement vector components)^ the predicted thickness 
errors are considerably larger than when the AT errors alternate 
in sign from one measurement vector component to another since 
the effects then tend to cancel . 

Calculations were also performed to test the performance of 
this data processing method for various aspect angles and both 
polarizations. Typical results are shown in Table 7. It was 
shown earlier in Figure 12, that the power reflection coefficient 
for vertical polarization at the Brewster Angle shows no dependence 
on the ice layer thickness so that the identification method fails 
completely. This is shown clearly in Table 7. 
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TABLE 7 

Minimum-Distance Pattern Classification 

5 K° Constant Temperature Measurement Error (8 Frequencies) 

Variable Angle of Incidence 

f = 1.80(.07)2.29 GHz 

Average Ice Thickness (Error in cm) 


Aspect Angle 

H.P. 

V.P. 

0° 

0.000 

0.000 

5° 

0.000 

0.000 

10° 

0.000 

0.000 

15° 

1.045 

1.201 

20° 

0.618 

1.101 

25° 

0.000 

0.000 

30° 

0.000 

0.000 

45° 

0.000 

0.910 

60° 

0.000 

43.543 

75° 

0.000 

0.003 


The data processing approach just described is a simple, brute- 
force approach which depends on accurate, absolute measurements 
of the apparent brightness temperature. Thus^ it is not well 
adapted to take into account small deviations in the apparent 
brightness temperature which may arise for example, from var- 
iations in the relative complex dielectric constant of the ice 
or variations in ice thickness. In the next section, we consider 
the construction of an alternate training set which will permit 
a minimum distance algorithm to be applied to ice models of 
unknown loss tangent, thickness distribution, etc. 


B. Fourier Series Fit 


It may be seen (from Figure 8, for example) that the reflection 
coefficient (which controls the brightness at the higher fre- 
quencies) is almost a periodic function of frequency for any given 
ice thickness. Thus the reflection coefficient, or the brightness 
temperature may be written as a Fourier series in frequency, e.g. 
at normal incidence 


Tjj(d) = Aq + J A^ cos n 27 TTf + J sin n2TTxf (5-15) 


where the “period" is l/twith T=2d/eT/c. (In fact, x is the 
two way propagation time of an impulseHhrough 
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the ice layer, and the A. are related to the reflection coefficients 
for the air-ice and ice-water boundaries. See Chapter II and the 
Appendix .) The brightness temperature may then be characterized 
by the parameter t and the Fourier coefficients A . If the 
brightness were measured at enough frequencies to Satisfy the 
sampling theorem for the thickest ice to be encountered, then the 
parameter could be unambiguously determined. If, however, fewer 
frequencies are utilized, then some ambiguities may exist in 
the estimate of x with corresponding errors in the estimated 
thickness. 

In the system considered here, the measurement vector is 
again a set of brighness temperatures on a span of frequencies 
The thickness estimate is made by assuming a value of thickness, 
and thus x, and using the known coefficients A to compute the 
corresponding brightness temperatures, so as t8 obtain an estimated 
of the measurement vector for the assumed thickness. The procedure 
is repeated for all thicknesses (1 cm - 100 cm) and the estimate 
thickness is taken as that for which the estimated vector is at 
a minimiin distance from the measured vector. The advantage of 
this procedure is that the effects of varying lossiness or of 
varying ice thickness can be accommodated by changing the Fourier 
coefficients. In comparison to the procedures of the next chapter 
the method requires considerably more computation. (That is, 
for each measurement, one might need up to 100 computations of 
brightness; in practice, one would utilize algorithms to restrict 
the number of trial thicknesses to a few in the neighborhood of 
the last thickness estimate, except for initial or ice boundary 
points). The memory requirements are much less, however, since 
only a few coefficinets need be stored. 

In the calculations reported here, the Fourier coefficients 
were determined by numerical harmonic analysis of curves such as 
those in Figure 9 at frequencies high enough that the galactic 
noise contribution could be considered negligible. Analyzed 
were curves representing aspect angles of 0°, 15^", and 30° for 
horizontal polarization and ice layer thicknesses of 10 cm, 

50 cm, and 90 cm. The results of the harmonic analysis are 
given in Table 8. The power reflection coefficient itself is 
approximately a cosine function, but since the apparent brightness 
temperature is proportional to (1 - |Rl ), the period of the 
temperature waveform to be analyzed was defined from one minimum 
to the next. At frequencies around 1 GHz or higher, the temper- 
ature curves are very much like true periodic functions, but not 
exactly, si nce^ though^ the galactic and atmospheric noise contribu- 
tion is very small compared to the physical temperature of ice, 
it is not zero. For this reason, the two minimum points on the 
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apparent brightness temperature curve for defining the period to 
be analyzed, could not be selected exactly, thus causing small 
errors in the computed harmonic amplitudes. The B's in Table 8 
are the amplitudes of the sine functions in the Fourier series 
representation and for lossless ice they all should be in 
principle zero. In Table 8 it can also be seen that the 
amplitudes of the “dc** term and the first harmonic tend to increase 
with increasing ice thickness. This is due to the fact that 
the ice model has small but non zero loss. Only if the ice 
layer were completely lossless, would the apparent brightness 
temperature harmonic amplitudes be independent of thickness. 

The coefficients for the three different thicknesses were 
averaged, and the "B were discarded to give the results of 

Table 9, 


TABLE 8 

Apparent Brightness Temperature Harmonic Analysis Results* 

(An. Bn in »K) 


o 

O 

It 

CD 

''o 

''l 

^2 

^3 

A4 

d=10 cm 

144.53 

-52.61 

10.18 

-1.91 

0.39 

d=50 cm 

145.33 

-52.86 

10.25 

-1.94 

0.39 

d=90 cm 

146.49 

-52.98 

10.10 

-1.92 

0.37 

0=15° 

d=10 cm 

142.96 

-54.29 

11.07 

-2.15 

0.47 

d=50 cm 

143.88 

-54.57 

11.08 

-2.19 

0.46 

d=90 cm 

146.50 

-54.29 

9.62 

-1.98 

0.24 

0=30° 

d=10 cm 

138.81 

-59.57 

13.64 

-3.04 

0.72 

d=50 cm 

140.31 

-59.78 

13.15 

-3.04 

0.62 

d=90 cm 

140.73 

-60.04 

13.66 

-3.06 

0.71 


*Based on the ice model of Chapter III, i.e. e. = 3.21 - j.0009/f. 
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TABLE 8 (continued) 


o 

O 

II 

^1 

^2 

^3 


d=10 cm 

1.71 

-0.31 

0.31 

0.08 

d=50 cm 

0.50 

-0.09 

0.09 

0.02 

d=90 cm 

0.56 

-0.16 

0.08 

0.00 

CD 

II 

CJ1 

o 

d=10 cm 

2.42 

-0.63 

0.41 

0.04 

d=50 cm 

0.44 

-0.78 

0.09 

0.02 

d=90 cm 

-5.56 

2.04 

-0.61 

0.13 

6=30° 

d=10 cm 

2.75 

-0.91 

0.53 

-0.01 

d=50 cm 

3.56 

-1.48 

0.57 

-0.12 

d=90 cm 

0.81 

-0.31 

0.14 

-0.02 


TABLE 9 

Averaged Apparent Brightness Temperature Fourier Series Coefficients 



^0 


^2 

^3 

A4 

9=0° 

145.45 

-52.82 

10.17 

-1.92 

0.38 

o 

LO 

II 

CD 

144.44 

-54.38 

10.59 

-2.11 

0.39 

0=30° 

139.95 

-59.80 

13.48 

-3.02 

0.68 


Then t is iterated over values corresponding to ice thickness from 
1 cm to 100 cm for each measurement vector, and the minimum-distance 
criterion of the previous sections is applied at each step to the 
measurement vector and to the training set vector produced by the 
Fourier series model. The computation time can be reduced, as 
mentioned above, by assuming that the iteration of x can be con- 
fined to some region bracketing the thickness of the preceding 
prediction. 

This method was applied to a number of different frequencies, 
three different scan angles, and different constant temperature 
errors (both positive and negative) up to 12 K°. The results are 
shown in Figures 35 (a) to (f) and 36 (a) to (f). For each scan 
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TEMPERATURE MEASUREMENT ERROR UK*} 



Figures 35 (a-b). 


Ice thickness prediction errors 
as a function of measurement 
errors, using five frequencies. 
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Figures 35 (c-d). Ice thickness prediction errors as 

a function of measurement errors, 
using five frequencies. 






TEMPERATURE MEASUREMENT ERROR (±K*) 



Figures 35 (e-f). Ice thickness prediction errors as 

a function of measurement errors, 
using five frequencies. 
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TEMPERATURE MEASUREMENT ERROR (iK*) 



TEMPERATURE MEASUREMENT ERROR (♦K") 


Figures 36 (a-b). Ice thickness prediction errors as 

a function of measurement errors, 
using six frequencies. 
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AVERAGE THICKNESS ERROR (cm) AVERAGE THICKNESS ERROR (cm) 



TEMPERATURE MEASUREMENT ERROR (±K«) 



TEMPERATURE MEASUREMENT ERROR (±K*) 


Figures 36 (c-d). Ice thickness prediction errors as 

a function of measurement errors, 
using six frequencies. 
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Figures 36 (e-f). Ice thickness prediction errors as 

a function of measurement errors, 
using six frequencies. 
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angle there are two curves. One corresponds to constant positive 
measurement errors, the other to constant negative measurement 
errors. The trends in this set of figures are similar to those 
of the preceding section. If the measurement errors are large 
{say 10K°, or more) then the average thickness errors can be 
reduced by using additional frequencies, as can be seen by 
comparing Figure 35 (a) and 36 (a), for example. Comparison of 
Figure 35 (a) with Figure 35 (f) shows that as a given frequency 
sample set is moved up in frequency, the ice thickness errors 
increase. If the error AT of Equation (5-14) is random with 
unifonn probability distribution between -10K° and +10K° instead 
of systematic, then the ice thickness errors are greatly reduced, 
as shown by the location of the X's in Figures 35 (a) and 36 (a). 

In general, the performance of this data processing method is 
quite similar to the one of the preceding section, as it should be, 
since the minimum criterion is the same. The advantage of the 
method is the ability to iterate over different ice types or thick- 
ness distributions by changing the value of the Fourier coefficients. 
For example, if the ice is multi-layered, the "dc" level remains 
about the same, but the amplitudes of the remaining harmonics are 
reduced. On the other hand, (see Figure 27 (b)), one observes that 
the increased loss due to inclusions in the ice causes the "dc** 
level to increase and the amplitudes of the harmonics to decrease. 

Since such variations in the apparent brightness temperature 
can easily be accounted for by changing the coefficients A ^the 
method can process data from ice with thickness variations"on the 
order of five centimeters, or from ice that has some loss; by iterating 
over the Fourier coefficients as well as x, the computation time, 
however, for the case of very lossy ice ('^^ 4 db/m) together with 
large multiple ice thicknesses ('^ ±. 5 cm) becomes excessive. 

The two high frequency data processing methods presented in 
this chapter depend heavily upon accurate, absolute apparent 
brightness temperature measurements. For this reason, they would 
not be suitable for measuring ice that has high loss and/or large 
multiple thicknesses, since both conditions cause the apparent 
brightness temperature to vary considerably. More sophisticated 
methods are needed to cope with the worst case possibilities in 
ice conditions. The next chapter presents two such methods. 
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CHAPTER VI 


DATA PROCESSING FOR LOSSY MULTI-THICKNESS ICE 
0. Introduction 


If the systems requirement is to measure smooth, clear ice, then 
the results of the previous chapter show that it is quite feasible 
to do so with five or six frequencies around 1 GHz or even slightly 
higher. If however, the measurements must be made over ice which 
might have inclusions or was of variable thickness, such as that 
commonly found on the Great Lakes, then the frequency range has to 
be chosen lower, and more complete data processing methods must be 
used. 


As to the choice of frequency span, consider for example Figures 
26 (a) and (b). If there were ± 5 cm multi-thickness present, the 
apparent brightness temperature would show no dependence on fre- 
quency or ice layer thickness at frequencies around 0.75 GHz. Thus, 
the frequency samples have to be chosen below 0.7 GHz. The upper 
frequency limit then is set by the multi-layering. The lower 
limit is set by the distortions in the apparent brightness temper- 
ature due tOpthe increasing galactic temperature contribution. 

For T = 2/r * , the lower limit appears to be about .3 GHz. The 
assumption of somewhat higher galactic temperatures would increase 
the lower frequency limit to 0.4 GHz or 0.5 GHz. This then (say 
0.3 GHz to 0.75 GHz) is the approximate frequency span available 

for measurement of lossy, variable thickness ice. 

As to the data processing system, the minimum-distance algorithms 
of the previous chapter were brute force methods that depended on 
a comparison of accurate absolute measurements of the apparent 
brightness temperature with a specific training set. Thus, they 
could not process temperature measurements from ice with losses or 
multi-layers. Unless such conditions were known in advance, the 
Fourier Series method could be adapted to handle unknown ice condi- 
tions but required much computation time since the training set 
vectors had to be recomputed for each measurement. This suggests 
that the Minimum-Distance method should have its training set 
expanded to range over all expected loss and variable thickness 
conditions. The memory requirement would be increased, but the 
computation time would be small, since no recomputation of the 
training set vectors would be necessary. 
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1. Mini mum- Pi stance Pattern Classification 

From the foregoing discussion, it is clear that to measure 
the thickness of lossy, multi-layered ice, the apparent brightness 
temperature measurements have to be made at frequencies lower than 
those considered in the previous chapter. In order to evaluate 
typical system performances at these lower frequencies, a number 
of the previous calculations were repeated, still assuming good 
ice of uniform thickness for the measurement vectors. The constant 
temperature error AT was taken to be 10°K. The minimum-distance 
algorithm was used as before. The results of these calculations 
are shown in Table 10. It can be seen that five temperature 
measurements at unequally spaced frequencies spanning 0.5 GHz to 
0.8 GHz will give acceptable results for constant temperature 
measurement errors of 10K°. It has been shown, however, that the 
variation in the apparent brightness temperature due to losses in 
the ice or multi-layers can be much greater than 10K°. 
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TABLE 10 

Average Thickness Errors For Five Frequency System 


Minimum Distance Pattern Classification 
(10°K Constant Temperature Error, Training Set And 
Measurement Set Based on Low Loss Ice of Constant Thickness) 


Frequency (GHz) 


Average ice 
Thickness 
Error (cm) 


.50 

.51 

.52 

.53 

.54 

.55 

.56 

.57 

.58 

.59 

.60 


X 

X 

X 

X 

X 







9.925 


X 

X 

X 

X 

X 






10.563 






X 

X 

X 

X 

X 


11.045 







X 

X 

X 

X 

X 

10.661 




X 

X 

X 

X 

X 




10.690 

X 


X 


X 


X 


X 



7.075 



X 


X 


X 


X 


X 

7.219 

X 


X 



X 



X 


X 

4.091 

X 

X 




X 




X 

X 

4.907 

X 




X 

X 

X 




X 

6.372 

X 



X 


X 


X 



X 

4.766 

.50 

.52 

.54 

.56 

.58 

.60 

.62 

.64 

.66 

.68 

.70 


X 

X 

X 

X 

X 







7.075 


X 

X 

X 

X 

X 






7.219 






X 

X 

X 

X 

X 


7.932 







X 

X 

X 

X 

X 

7.085 




X 

X 

X 

X 

X 




7.543 

X 


X 


X 


X 


X 



1.236 



X 


X 


X 


X 


X 

1.035 

X 


X 



X 



X 


X 

1.050 

X 

X 




X 




X 

X 

2.354 

X 




X 

X 

X 




X 

2.166 

X 



X 


X 


X 



X 

1.420 

.50 

.53 

,56 

.59 

.62 

.65 

.68 

.71 

.74 

.77 

.80 


X 

X 

X 

X 

X 







2.317 


X 

X 

X 

X 

X 






2.962 






X 

X 

X 

X 

X 


2.998 







X 

X 

X 

X 

X 

3.322 




X 

X 

X 

X 

X 




2.606 

X 


X 


X 


X 


X 



1.219 



X 


X 


X 


X 


X 

1.216 

X 


X 



X 



X 


X 

1.015 

X 

X 




X 




X 

X 

1.372 

X 




X 

X 

X 




X 

1.196 

X 



X 


X 


X 



X 

0.779 
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TABLE 10 (continued) 
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A more realistic simulation of system performance, in which 
the training vectors were based on good ice of constant thickness, 
and the "experimental" results were based on good ice of variable 
thickness, is shown in Table 11. For this case, six frequencies 
were required for reasonably good results. The case of a constant 
10K° measurement error and ice of uniform thickness has been 
included for comparison purposes. It can be seen that the addition 
of the sixth frequency helps to reduce the average error in 
predicted ice thickness. Note that the thickness variations 
have been assumed to lie on top of the basic ice layer in this 
set of calculations. Thus, if a measurement vector is constructed 
from the apparent brightness temperature for a 70 cm ice layer 
with 5 cm multi-layers on top, the predicted thickness should be 
the average of all the layers, i.e. 72.5 cm, or 2.5 cm more than 
the basic ice layer thickness. (This can be seen in the average 
thickness predictions in Table 11 when an average ice thickness 
prediction "error" of 0.5 cm is expected for the 1.0 cm multi- 
layering, etc..) It is clear that six frequencies can be chosen 
to give good performance for up to 5 cm of multi -layering, and 
it is suggestive for system design that best results are usually 
obtained when the total span of frequencies is as large as 
possible with the frequencies irregularly (rather than uniformly) 
spaced. The method fails (with six frequencies) for 10 cm 
multi -layering. If in addition, the effects of realistic loss 
tangents were incorporated into the measurement vectors, the 
performance would deteriorate even more. 

The Minimum-Distance method of data processing can be extended 
to accommodate some of the problems caused by ice of variable 
loss or thickness by enlarging the training set in the memory. 

This training set enlargement has been shown schematically in 
Figure 37, which depicts an array of 9 training sets for three 
different loss tangents and three different thickness variations 
(± 1 cm to ± 4 cm). The training set now has 900 six-dimensional 
vectors, assuming a six frequency system. (In order to minimize 
storage requirements for the computations to follow, each 
training set is based on 100 rather than 199 thicknesses in 
the interval 0 to 100 cm. Clear ice is modeled by e. = 3.20 - 
j.OOl/f, where f is in GHz. More lossy ice, such asHhat measured 
by Vickers (1975)^ is represented by e. = 3.10 - j.04/f and e. = 

3.00 - j.08/f, corresponding to attenlation of about 1 db/m Jnd 
4 db/m respectively.) 

If the measurements were taken from one of the nine types 
of ice in the training set without temperature error, then 
a perfect fit could be found for each measurement vector. A 
realistic test of this data processing procedure would be to 
assume that the ice to be measured was different from any of the 
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TABLE n 

Average Thickness Errors For Six Frequency System 

Minimum Distance Pattern Classification 
(Training Set Based on Low Loss Ice-Constant Thickness. 
Measurement Set Based on Low Loss Ice-Variable Thickness.) 


Frequency (GHz) 


Average Ice Thickness Error (cm) 


.50 .51 


.52 .53 .54 .55 .56 .57 


.58 .59 .60 


X 


X 


o 

<y> 


X 


X 

X X 


X X 


XXX 

XXX 


X X 

X X 


X 


X 


X 


X 


X 


X 


X 

X 

X 

X 

X 


.50 .52 .54 .56 .58 .60 


.62 .64 .66 .68 .70 


X 



XXX 
X X 


X X 


X 


X 


X 

X 


X 

X 


X 

X 

X 

X 


X 

X 

X 


X 

X 


X 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

X 


10 K° 1 cm 

9.621 0.495 
9.646 0.495 
8.872 0.495 
6.553 0.495 
4.819 0.495 
4.613 0.495 
4.704 0.495 


4.613 0.495 
4.560 0.495 
5.025 0.495 
0.930 0.495 
0.606 0.495 
0.422 0.495 
0.558 0.495 


5 cm 10 cm 

10.530 30.294 
11.859 30.749 
12.568 32.239 
5.508 32.367 
4.131 25.997 
3.766 35.231 
4.168 34.869 


3.766 35.231 
4.794 35.970 
5.905 36.525 
4.178 30.033 
2.731 34.442 
2.601 33.239 
2.550 36.626 


X 


X 


TABLE n (continued) 


Frequency (GHz) Average Ice Thickness Error (cm) 


.50 

.53 

.56 

.59 

.62 

.65 

.68 

,71 

.74 

.77 

.80 

10 r 

1 cm 

5 cm 

10 cm 

X 

X 

X 

X 

X 

X 






1.178 

0.495 

2.967 

35.460 



X 

X 

X 

X 

X 

X 




1.294 

0.495 

3.691 

37.166 






X 

X 

X 

X 

X 

X 

1.608 

0.495 

4.764 

35.153 

X 

X 

X 






X 

X 

X 

1.146 

0.495 

4.977 

37.457 

X 

X 



X 

X 




X 

X 

0.862 

0.495 

3.033 

37.854 

X 


X 


X 


X 


X 


X 

0.812 

0.495 

6.678 

23.038 

X 

X 

X 



X 



X 


X 

0.101 

0.495 

2.950 

36.412 

.50 

.54 

.58 

.62 

.66 

.70 

.74 

.78 

.82 

.86 

.90 





X 

X 

X 

X 

X 

X 






0.422 

0.495 

2.601 

33.237 



X 

X 

X 

X 

X 

X 




0.495 

0.495 

3.626 

37.894 






X 

X 

X 

X 

X 

X 

0.525 

0.495 

6.113 

38.279 

X 

X 

X 






X 

X 

X 

0.407 

0.495 

4.739 

36.945 

X 

X 



X 

X 




X 

X 

0.106 

0.495 

4.420 

32.168 

X 


X 


X 


X 


X 


X 

7.975 

0.495 

13.143 

22.010 

X 

X 

X 



X 



X 


X 

0.131 

0.495 

3.623 

21.839 

.50 

.55 

.60 

.65 

.70 

.75 

.80 

.85 

.90 

.95 

1.00 





X 

X 

X 

X 

X 

X 






2.259 

0.505 

4.000 

32.075 



X 

X 

X 

X 

X 

X 




1.947 

0.505 

5.560 

35.314 






X 

X 

X 

X 

X 

X 

2.068 

0.575 

9.465 

35.804 

X 

X 

X 






X 

X 

X 

2.274 

0.756 

8.405 

23.063 

X 

X 



X 

X 




X 

X 

2.641 

0.505 

9.827 

28.412 

X 


X 


X 


X 


X 


X 

19.892 

0.932 

23.796 

26.942 

X 

X 

X 



X 



X 


X 

3.437 

0.585 

4.286 

32.580 





3.20- j.OOI/f 



3.l0-J.04/f 


3.00- j.08/f 


Figure 37. Enlarged training set used for ice of variable loss 
and thickness. 




nine types in the memory and the measurements were then corrupted by 
some temperature error AT. 

In the example considered here, the measurement vectors are 
assumed to come from ice represented by e. = 3.05 - j.06/f with 
a variable thickness of ± 2 cm. In FigurJ 37, this type of ice 
would be located by the matrix I,. When arbitrary constant temper- 
ature errors were added, it was found that errors appeared in the 
predicted ice thickness in a random fashion depending on the 
size of the constant temperature error. These random errors in 
the predicted ice thickness which occur when there is a systematic 
measurement error^ can be minimized by "randowiitij^g*' them. .The 
procedure is to add six different temperature errors (± 10K°, + 15K°, 
+ 20K°) to each measurement vector, make the six corresponding ice 
thickness predictions, and then take a "vote" among the six different 
predicted ice thicknesses. By choosing the prediction closest to 
the average of the six predictions , this rather 

than smoothes out, spurious errors, and a typical result is shown 
in Figure 38. 



Figure 38. Ice thickness errors using minimum-distance algorithm 
with expanded training set. 
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This method of data processing obviously requires more storage 
and computation time but has the advantage of providing some in- 
formation about the lossiness of the ice and large scale thickness 
variations. 

The existence of constant (i.e. systematic) temperature errors 
of more than ± 20K° produces a significant deterioration in system 
performance. The reason, again, is that one is attempting to 
match absolute brightness chosen from a cosinusoidal variation of 
brightness with thickness in which the average ("dc") level and the 
amplitude depend on the thickness variations and the losses in the 
ice. It would thus be desirable to develop a method of data pro- 
cessing which would depend on the relative brightness in the several 
frequency channels rather than the absolute brightness temperature. 


2. ' Temperature Quantization 

The minimum-distance method of the preceding section can handle 
systematic temperature errors on the order of + 15K°, even in the 
presence of brightness due to changes in the ice properties, but 
fails for larger errors. To develop a data processing method that 
is independent of constant measurement errors, one may exploit the 
fact that the periodicities, and thus the relative levels, of the 
brightness vs frequency curves are not much affected by losses or 
multiple thickness effects. 

To implement such a system, one must again measure the 
brightness at a number of different frequencies. At each pair 
of frequencies, the ice will characterize the fact that the bright- 
ness at one frequency is larger, smaller, or equal to the brightness 
at the other frequency. The result of this decision, or quantiza- 
tion, may now be represented by a ternary number (e.g. "less than" 

= 0; "greater than" = 1; "equal to" =2). If n different frequencies 
are utilized, n(n-l)/2 such comparisons may be made, and the 
"measurement vector" will thus be a sequence of n(n-l)/2 ternary 
numbers. The training sets (one set for each type of ice, i.e. 

9 sets for example of Figure 37) will each consist of, say, 100 
such sequences, one for each possible thickness between 1 cm and 
100 cm. The criterion for a best fit between the measurement 
vector and a member of the training set is based on the number of 
matches (rather than some minimum distance) in the two ternary 
sequences. 

In practice, two modifications are needed to successfully 
implement this scheme. In the first place, to accommodate random 
errors and to maximize the information content of the sequence, 
the criterion for "equality" must be modifed to include a "dead- 
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band" z\T such that temperatures T, , J 2 are considered "equal" if 
1^1 ’ ^ Clearly, as at' is increased, the number of 

2‘i in the characterizing sequence^il 1 increase from almost none 
to almost all. Secondly, it was found that, in order to minimize 
the effect of adventitious matches, it was desirable to impose a 
more stringent criterion for a match. In the algorithm adopted 
here, a match is considered valid only if the thickness derived 
from the two best matches are within 2 cm of each other. 

In the actual implementation of the procedure, seven 
frequencies, representing 21 comparisons, were found to give 
good results. The search for a valid thickness estimate is 
begun by assuming a particular AT (in the example of Table 11, 
the assumed value was aT = 10°K). The measurement vector based 
on this aT^ is then compared with the training set (based on 
the same aT ) for one type of ice. If a "valid" match is 
found, the thickness is taken to be the average thickness of the 
two best matches. If a valid match is not found, the procedure 
is repeated for AT = 5°K, then (if necessary) for AT = 15®K, 
and finally (if necessary) for AT = 2°K. This proced0re either 
yields a thickness prediction or it does not. 

The entire process is repeated for each of the 9 ice types, 
so that one obtains up to 9 thickness estimates. The final 
estimate of thickness is made by taking a "vote" i.e. by 
choosing the thickness most often represented among the nine 
results . 

To illustrate the effectiveness of this system, it was tested 
for the ice type represented in Figure 37 by matrix K (e = 3.05 
- j.06/f, with thickness variation of ± 3.5 cm). The^results are 
illustrated in Figure 39, which shows the average thickness error. 
As can be seen from the figure, the random appearance of incorrect 
thickness prediction has been practically eliminated by the use 
of the voting procedure. 

Although the data storage requirement is rather large 
(9 X 4800 21 digit ternary numbers for 4 levels of AT and 
a 3x3 ice type array), the data processing consists onfy of simple 
digit by digit comparisons. Because systematic errors (introduced, 
for example, by changes in the cosmic background or by calibration 
drift in the equipment) do not effect this temperature quantization 
procedure, it appears to be the most useful of those considered. 


Ill 




Figure 39. Ice thickness prediction errors using temperature 
quantization method with expanded training set. 
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CHAPTER VII 


SEA ICE THICKNESS MEASUREMENT 


0. Introduction 


There is now considerable interest in the ability to measure 
the thickness of sea ice. In part, this interest has been gener- 
ated by results of recent climate models which indicate that the 
distribution and the amount of sea ice has an important influence 
on world climate, and in part, because such information is necessary 
for successful surface and submarine navigation in the polar areas. 

Many different approaches have been taken towards solving the sea ice 
thickness measurement problem. They include radar scatterometry, 
Parashar (1974), side-looking airborne radar. Rouse (1969), radar 
utilizing holographic principles, lizuka (1971), satellite IR data 
analysis, LeSchack (1975), inference of thickness distribution from 
large scale geophysical properties, Thorndike (1975), and measure- 
ment of microwave emissions from aircraft and satellites. Strong 
(1970), Adey (1972), Gloersen (1972), Wilheit (1972), Campbell (1974), 
Gloersen (1974), Meeks (1974), Tuiri (1974), Gloersen (1975), Tiuri 
(1975), Some of these methods and a laser profilometer were combined 
in a comprehensive measurement program during the passage of the 
oil tanker Manhattan, Biache (1971). As for qualitative systems, 
it should be noted that radiometric methods have been used successfully 
for some time to detect icebergs in the ocean, Ladely (1967), Sperry 
(1967). A number of both active and passive microwave systems have 
successfully distinguished between different ice categories (first 
year, multi-year) without being able to provide quantitative thickness 
data. The radar system described by lizuka (1971) can measure sea 
ice thickness up to 4 cm with 9% accuracy. The radiometric systems 
at 765 MHz of Adey (1972) and at 600 MHz of Tiuri (1974), (1975) 
have shown a direct, but qualitative relationship between the 
apparent brightness temperature and sea ice thickness. It is now 
clear that the limited quantitative success of most of these systems 
is due primarily to the complex nature of salt water ice in the 
natural state. 

Since we have shown, in the previous chapters, that it is 
feasible not only to measure but to map the thickness of fresh 
water ice layers using radiometric methods, it is natural to 
applying the same methods to sea ice. Because the electrical 
properties of sea ice are significantly different from those of 
fresh water ice, we first review the extensive literature on this 
subject, before estimating system performance. 
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1. Literature Review 


The primary difference between salt water ice and fresh water 
ice is the presence of brine which has been trapped in the salt 
water ice by the freezing pockets. The size, number, and 
orientation of these pockets are the major source of variability 
in the properties of sea ice. Furthermore, the brine pockets 
are not stationary, Kingery (1963). Brine drainage occurs under 
the influence of gravity through interconnected channels in the ice, 
and the brine pockets themselves diffuse through the ice toward 
the warm side (i.e., the bottom of the sea ice layer) under the 
influence of a temperature gradient. Thus, through these two 
different mechanisms, the salinity profile of a layer of sea ice 
will change considerably over a period of time. It is this 
presence, or absence, of brine in sea ice which makes it possible 
to distinguish new ice from old ice radiometrically. The higher 
brine content in the new ice tends to make it more lossy, hence 
more emissive, thus giving it a higher apparent brightness temper- 
ature in the microwave frequencies. An extensive discussion of 
the radiometric properties of sea ice at X-Band and above has been 
given by Poe (1972). 

The amount of brine entrapped by the freezing process is a 
function of the rate of freezing. Sea water suddenly exposed to 
temperatues of -40°C may entrap enough brine to give a 2% 
salinity, Pounder (1965). By way of comparison, sea water has 
3.5% salinity. Usual ly, young ice has surface salinities of 
,8% to 1.0%, Bennington (1963). The orientation of the brine 
inclusions with respect to the incident field is also very 
important. It has been found that the loss factor changes by 
as much as a factor of four, depending on the relative orientations 
of the brine inclusions and the electric field, Sackinger 
(1972). Sea ice anisotropy has also been measured by Campbell 
(1974) using a short pulse radar system. Thus, the electrical 
properties of salt water ice are a function of how fast the ice was 
frozen, how long it has been frozen, what temperature gradients 
it has experienced, and the physical orientation of the brine 
pockets with respect to the electric field. An additional 
variable is the physical temperature of ice at the time of 
measurement. This, in part, may explain the great variation 
in reported loss measurements of salt water ice. Furthermore, 
most of the results to date come from laboratory measurements 
made on artificial sea ice water, and it is possible that the 
laboratory-grown ice does not represent accurately sea ice in 
the natural state because of slightly different brine distribu- 
ti ons . 
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To illustrate the great range of loss values possible, some 
of the results reported in the literature. Cook (I960), Iglesias 
(1967), Addison (1969), Gudmandsen (1971), Hoekstra (1971), lizuka 
(1971) are listed in Table 12. 


TABLE 12 

SALT WATER ICE PROPERTIES 


f(GHz) 

db/m 

Ref 

°C 

Salinity % 

.10 

26. 

Addison (1969) 

-12.5 

0.512 

.10 

4.1 

Addison (1969) 

-33. 

0.512 

.10 

23. 

Addison (1969) 

-15. 

1.25 

.10 

7.7 

Addison (1969) 

-35. 

1.25 

.10 

20. 

Addison (1969) 

-15. 

1.58 

.10 

7.1 

Addison (1969) 

-35. 

1.58 

.10 

23. 

Cook (1960) 

-10. 

0.70 

.10 

0.75 

Cook (I960) 

-40. 

0.70 

.10 

10. 

Cook (I960) 

-10. 

0.35 

.10 

0.77 

Cook (I960) 

-40. 

0.35 

.45 

.62 

Gudmandsen (1971) 

-20. 

0.2 

.45 

55.0 

Gudmandsen (1971) 

-20. 

0.4 

.50 

108. 

Hoekstra (1971) 

-10. 

0.8 

1.00 

135. 

Hoekstra (1971) 

-10. 

0.8 

1.00 

2.01 

Iglesias (1967) 

-1. 

0.002 

1.00 

0.06 

Iglesias (1967) 

-60. 

0.002 

2.00 

129. 

Hoekstra (1971) 

-10. 

0.8 

10.00 

120. 

Gudmandsen (1971) 

-20. 

0.4 


These losses are reported in many forms, e.g. nepers /m, db/m, loss 
tangent, imaginary part of'the complex dielectric constant; so 
for ease of comparison, all the results have been converted to 
db/m. In some cases^ an assumption had to be made regarding the 
real part of the complex dielectric constant. It is clear from the 
table that there is a very large variation in the reported results, 
but that in general the losses are significantly higher than those 
found for fresh water ice. 
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Another complication in the salt water/salt ice system is 
the nature of the boundary between the water and the ice layer. 

In fresh water this boundary usually is well defined. In labor- 
atory experiments on salt water freezing, it has been found 
that icicles grow downward from the ice layer into the water, 
Palosuo (1963). These icicles tend to form an impedance matching 
layer between the ice and the water. Thus, the oscillations in 
the apparent brightness temperature due to interference between 
the two boundaries would be further reduced. Other complications, 
such as small and large scale roughness on the top surface and 
air inclusions and cracks in the ice layer make quantitative 
interpretation of apparent brightness temperatures of sea ice 
even more difficult. 


2. ^parent Brightness Temperature Calculations 

Because of the large attenuation likely to be encountered in 
sea ice, one would expect that the ripples in the apparent brightness 
temperature would be highly damped, and the brightness temperature 
would increase assymptotically to some value given by the product 
of the ice emissivity and the physical temperature of the ice. A 
number of calculations have been made using the apparent brightness 
model developed previously with representative values of attenuation 
for the ice and assuming a salinity of 35 ppt for the underlying 
sea water. ^ 


Note that^for the frequency ranges considered here (0.1 GHz 
to 2.0 GHz)^ the complex dielectric constant of fresh water does 
not change much, whereas, for typical sea water (salinity^ 35 ppt). 
the imaginary part of the complex dielectric constant increases 
rapidly at low frequencies due to the ionic conductivity of the 
dissolved salts. (See Figures 40 (a) and 40 (b)). Representative 
values of the dielectric constant and the power reflection coefficient 
of fresh and salt water are given in Table 13. 


TABLE 13 


REFLECTION COEFFICIENT OF FRESH AND SALT WATER 


f=0.1GHz 

f=1.0GHz 


s=0 ppt 

e 

W 



87.7- J0.9 0.651 

86.7- J9.1 0.651 


e 


w 


s=35 ppt 


R 


2 


75.6-J523.2 0.877 

74.8-J59.7 0.679 
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FREQUENCY (GHx) 


Figures 40 (a-b). Relative complex dielectric constant of fresh 

and salt water. 
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The power reflection coefficient as a function of ice thickness, 
with ice attenuation as parameter, is shown in Figures 41 (a), (b), 
(c), (d) , and (e). The corresponding apparent briqhtness temperatures 
are shown in Figures 42 (a), (b), (c), (d), and (e). From this 
set of figures^ it can be seen that if the ice has an attenuation 
of approximately 8 db/m or less, the temperature quantization 
method of data processing could be utilized to measure Ice thickness 
up to a meter. If the attenuation is more than 32 db/m for example, 
only ice layers less than 10 cm to 20 cm thick could be success- 
fully measured by any data processing method that relied on estimates 
of the period of oscillation of the brightness. In practice, the 
tolerable attenuations are even less, since actual sea ice would 
probably exhibit considerable surface texture, some variability 
in thickness, and an ill-defined transition layer between the 
ice and the water. All three of these conditions tend to smooth 
out the ripples in the apparent brightness temperature. 



Figure 41 (a) . 


Power reflection coefficient as a function of ice 
thickness and ice loss. Salinity of water = 35 ppt. 
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Figures 41 (b-c). Power reflection coefficient as a function of ice 

thickness and ice loss. Salinity of water = 35 ppt. 
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Figures 41 (d-e). Power reflection coefficient as a function of ice 

thickness and ice loss. Salinity of water - 35 ppt. 
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Figures 42 (a-b). Apparent brightness temperature as a function of 

ice thickness and ice loss. Salinity of water 
= 35 ppt. 
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Figures 42 (c-d). Apparent brightness temperature as a function of 

ice thickness and ice loss. Salinity of water 
= 35 ppt. 
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ICE THICKNESS (cm) 


Figure 42 (e). Apparent brightness temperature as a function of 
ice thickness and ice loss. Salinity of water 
= 35 ppt. 


The change in the power reflection coefficient (as a function 
of frequency) with increasing attenuation is shown in Figures 43 
(a) and (b) for a 10 cm ice layer and in Figures 44 (a), (b), 
and (c) for a 100 cm ice layer. The apparent brightness temper- 
atures corresponding to these power reflection coefficients are 
shown in Figures 45 (a) and (b) and 46 (a), (b), and (c). This 
set of figures offers further evidence that no more than 8 db/m 
attenuation can be tolerated if ice thicknesses up to a meter 
are to be measured by the methods considered here. 

For very high loss, the apparent brightness temperature 
changes almost monotonically from that for sea water to that for 
thick ice. Thus, if the attenuation were known to be high, a 
measurement of brightness temperature at one or two frequencies 
would provide a fairly good estimate of thickness directly. 
However, once the total attenuation through the layer (i.e. the 
product of attenuation and thickness) reaches 15 db or so, 
the thickness has almost no influence on the reflection coeffic- 
ient. Thus, the brightness which is then determined by the 
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dielectric properties and physical temperature of the ice, not 
its thickness. Thus, the single frequency systems (Tiuri (1975), 
Adey (1972)) in the 600 MHz - 800 MHz show the expected relation 
between increasing brightness and increasing thickness, but can 
only provide qualitative estimates of the actual thickness. 
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POWER REFLECTION COEFFICIENT 


d « 10 cm 




Figures 43 (a-b). Power reflection coefficient as a function of 

frequency and ice loss. Ice Thickness = 10 cm, 
water salinity = 35 ppt. 
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Figures 44 (a-b). Power reflection coefficient as a function of 

frequency and ice loss. Ice thickness « 100 cm, 
water salinity = 35 ppt. 
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Figure 44 (c). Power reflection coefficient as a function of 

frequency and ice loss. Ice thickness = 100 cm, 
water salinity = 35 ppt. 
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Figure 45 (a). Apparent brightness temperature as a function of 
frequency and ice loss. Ice thickness = 10 cm, 
water salinity = 35 ppt. 
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Figure 45 (b). Apparent brightness temperature as a function of 
frequency and ice loss. Ice thickness = 10 cm, 
water salinity == 35 ppt. 



FREQUENCY (GHz) 


Figure 46 (a). Apparent brightness temperature as a function of 
frequency and ice loss. Ice thickness = 100 cm, 
water salinity = 35 ppt. 
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Figures 46 (b-c). Apparent brightness temperature as a function of 

frequency and ice loss. Ice thickness = 100 cm, 
water salinity = 35 ppt. 


129 


CHAPTER VIII 


SYSTEM CONSIDERATIONS 


0. Introduction 

The preceding results show that it is feasible to measure 
fresh water ice thickness if brightness measurements at a sufficient 
number of frequencies are available. It remains to investigate 
the practicality of such a system from the standpoint of the 
required receiver sensitivity, as determined by the integration 
times, bandwidths, and receiver noise figures. For a mapping sys- 
tem, one must also consider the effects of aircraft height and 
velocity as they affect the scanning process. It will be shown 
that a range of multichannel radiometer configurations within the 
current state-of-the art in radiometer hardware will satisfy the 
performance requirements for a feasible mapping system. 

1 . Bandwidth 

All the data processing schemes considered thus far have 
assumed a zero bandwidth. A radiometer, however, receives power 
proportional to the bandwidth, and therefore the effect of a 
finite bandwidth must be taken into consideration. Figure 47. 
shows the apparent brightness temperature as seen by radiometers 
at 1 GHz center frequency but with varying bandwidths. This 
calculation was performed by averaging the apparent brightness 
temperatures of a number of frequencies distributed about a 
center frequency, in this case 1 GHz. Thus the calculation 
simulates the effect of an ideal bandpass filter. The effect of 
other types of bandpass filter shapes is considered by Popov 
(1976). 

As the bandwidth is further increased, the apparent bright- 
ness temperature becomes modulated by the function sin (d Ag)/ 

(d A 3 ). ThuSj the effect of finite bandwidth is similar to that 
due to multiple thickness but poses much less of a problem in 
practice, because the bandwidths are small. It can be seen 
from Figure 47. that a 10 MHz bandwidth, for example, would have 
a negligible effect on the apparent brightness temperature. 

If larger bandwidths are required to improve the sensitivity 
or reduce the integration time, then the smoothing effect has 
to be allowed for in the data processing scheme. 

2. Receiver Sensitivity for Scanning Systems 

The sensitivity of a radiometer is closely related to the 
bandwidth. A well known expression, Kraus (1966), for the minimum 
detectable temperature AT is 
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Figure 47. The effect of finite bandwidth on the apparent brightness 
temperature. 



(8-1) 


AT = 


/aT 


where 


T = 290. (F - 1) + T . 
sys ' ' ant 


K is the system constant 
F is the receiver noise figure 
Tant antenna temperature 

Af is the receiver bandwidth 
T is the integration time. 


( 8 - 2 ) 


The system constant K depends on the nature of the system. 

A Dicke receiver with square wave modulation and a narrow band 
video amplifier, has K - 2.22, for example; a total -power re- 
ceiver has K = 1.0, Kraus (1966). A representative plot of the 
sensitivity as a function of the integration time for a total - 
power receiver is given in Figure 48. A noise figure of 4 and 
antenna temperatures of 100° K and 200° K have been assumed in 
the figure. In practice, Dicke receivers are usually employed 
since they greatly reduce the gain instabilities inherent in the 
total -power receiver, but even if the T of Figure 48. is increased 
by a factor of two, the sensitivity is still quite adequate in 
terms of the kinds of temperature errors (5°K - 10°K) which have 
been assumed in the thickness prediction studies of Chapter V and 
VI. 


In a practical mapping system, the available integration 
times are determined by the flight velocity and system geometry. 
In a thickness profiling (not scanning) mode, the integration 
time is just the time required to fly over the antenna footprint 
on the ground. A brief consideration of typical aircraft velocity 
and altitudes shows that for any practical beamwidth, the 
available integration time is adequate to meet a sensitivity 
requirement of a few °K. The implementation of a scanning system 
presents a design conflict however. Good ground resolution 
requires that the antenna beamwidth be small, which in turn 
requires that the scan rate be high. High scan rates imply 
short integration times causing the minimum detectable temper- 
ature AT to increase. Consider a system in which the aircraft 
flies at an altitude h, with velocity v, antenna beamwidth 0n, 
and a scan angle 0 (Figure 49.). If the axis of rotation is° 
tilted degrees from the vertical, it is easy to show that 
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Figure 48, Receiver sensitivity as a function of bandwidth and 
integration time. 



Figure 49. Scanning geometry. 
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(8-3) 


K T 






AT 


h Af 


This is an expression for the antenna beamwidth compatible with 
the specified receiver sensitivity AT. This beamwidth has been 
calculated as a function of velocity/height ratio and plotted 
in Figure 50. for representative values of noise figure, band- 
width, and receiver sensitivity. The calculations have assumed 
K = 2.22, 0 = 60°, (j) . = 90°, and T = 200°K. As an example of 
a typical flight system^ consider an^aircraft flying 100 mph, at 
2000 feet with a v/h ratio of 0.073 sec” . If the system noise 
figure is 4 and bandwidth is 10 MHz, then 0^ = 10° for AT = 1°K 
and 0g = 2° for AT = 5°K. Since the results^of the previous 



Figure 50. Required beamwidth as a function of velocity /height 
ratio. Curves are calculated assuming K = 2.22, 
e. = 80°, (j) . = 90°, and T, = 200°K. 

S u a 
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sections indicate that the proposed data processing systems can 
operate effectively with errors of the order of 5°K, it appears 
that for many systems of interest a beam scanning design can be 
implemented. For more severe opierational requirements (e.g. very 
high speed platforms or significantly narrower beamwidths), it 
may be necessary to use multiple beam systems to achieve the 
desired sensitivity. 

3. Receiver Configuration 

From the discussion of data processing methods^ it is clear 
that five, six, seven, or more frequencies are needed to predict 
the ice thickness. Furthermore, since the radiometer will be 
airborne, it is important to keep the integration time for each 
thickness measurement to a minimum. These two requirements 
dictate a parallel integration configuration, one example of 
which is shown in Figure 51. in which all channels are simultan- 
eously in operation. The choice of actual system implementation 
(e. g. whether to achieve channel separation in the r.f. section 
by using multiple local oscillators or by using i.f. amplifiers 
of different frequencies) will depend primarily on the availability 
of hardware, image rejection considerations, etc.. 
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Figure 51. One possible configuration of a 







IX. CONCLUSIONS 


This study presents the basic calculations from which the 
feasibility of a microwave radiometer system for mapping the thickness 
of a lake ice may be evaluated. The study is based on the well 
established principle that the reflection coefficient of an ice layer 
over water depends upon the thickness of the ice and the frequency of 
observation. The major emphasis in this study has been the development 
of mapping systems for the quantitative determination of ice thickness 
under realistic conditions. 

An analytical model for the apparent brightness temperature ns 
a function of ice thickness has been developed and elaborated to 
include such variables as galactic and atmospheric noise temperatures, 
aspect angle, polarization, temperature gradients in the ice, 
and the presence of transition layers above and below the ice. Look 
angles between 0° (normal incidence) and about 55° were found to 
give satisfactory performances; circular or 45° linear polariza- 
tion gave the best results at the larger look angles. The effect 
of temperature gradient on the apparent brightness temperature was 
found to be negligible. A snow layer on top of the ice was found 
to act as a matching layer or an absorber depending primarily on 
the moisture content of the snow; a few centimeters of dry snow 
could be easily tolerated for most of the systems considered here. 

A slush or a water layer on top of the ice changes the reflection 
coefficient radically. Water layers more than a few mm thick were 
found to degrade severely the operation of any measurement system. 

Small scale first surface roughness was found to change the apparent 
brightness temperature by only a few degrees, and thus^ would present 
no problem to thickness measuring systems. The presence of a range 
of ice layer thicknesses within the antenna footprint was found to 
place an upper limit on the usable frequency. A 10 cm variation 
in ice thickness, for example, limited usable frequencies to below 
0.7 GHz. The presence of air bubbles in the ice layer, which tend 
to make the ice more lossy, was found not to degrade the thickness 
estimates, at least for the range of attenuation typical of fresh 
water ice. 

The second objective of the study was the development of 
effective data processing methods. Two methods were developed 
to predict the ice layer thickness from brightness temperature 
measurements: the minimum-distance fit, and ternary temperature 

quantization. It has been found that the thickness of ice layers 
up to a meter can be estimated with reasonable accuracy using as 
few as three frequencies, but they must lie below 100 MHz, and the 
data processing methods must take into account the variations in 
galactic noise temperature. In addition, the design of an airborne 
antenna at such low frequencies is a formidable problem. At 
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higher frequencies (1.0 GHz or more)^ the antenna design problem 
is much more tractable, but as many as six unevenly spaced fre- 
quencies, spanning approximately .5 GHz, are needed to obtain 
accurate estimates of the thickness of clear uniform ice up to 
a meter thick. If the ice becomes lossy due to air bubbles 
and is of non-uniform thickness, then seven frequencies in the 
span .4 GHz - .7 GHz are required for the minimum-distance method 
The temperature quantization method provides performance 
comparable to the minimum-distance method but has the advantage 
that it is less sensitive to systematic temperature errors. 

A brief system study has been performed to evaluate the 
required receiver design for an airborne scanning system. The 
effects of antenna beamwidth, scanning rate, bandwidth, system 
noise figure^ and integration time were considered for typical 
vehicle velocity/altitude ratios. No problems are anticipated 
in meeting the receiver sensitivity requirements with state-of- 
the-art components. ThuS; it has been shown that a multichannel 
microwave radiometer can be used to map the thickness of lake ice 
subject to certain reasonable constraints upon the thickness 
variations and losses in the ice layer. 

The feasibility of using a similar system to estimate the 
thickness of sea ice was also investigated. Reported loss values 
for sea ice are large and highly variable, depending on the state 
and history of the ice. Calculations show that the apparent 
brightness temperature is a function of sea ice layer thickness 
only for very thin layers ( a few tens of centimeters). Radio- 
metric methods^ thus, can distinguish new ice from old ice since 
the brine in the old ice has partially drained out, but do not 
appear to be feasible for quantitative measurements of thickness 
of the order of a meter, unless the ice to be measured is of 
exceptionally low loss. 
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APPENDIX 


Brightness Temperature From Impulse Response 

It is well known from Fourier transform theory that a periodic 
spectral response corresponds to a set of impulse functions in the 
time domain and vice versa. Consider for example, the Fourier 
transform of two impulses in the time domain of arbitrary relative 
magnitudes, one at t=0, the other delayed in time by some constant 
T , i . e . 

E(jw) = ,^[s(t) + A6(t-x)] = 1.0 + Ae . (A-1) 

The square of the magnitude of the function in the frequency domain 
is 

iF(ju)) = 1.0 + + 2A cosuT . (A-2) 

The argument of the cosine function can be rewritten as 

|F{ju) 1^ = 1.0 + A^ + 2A cos ^ , (A-3) 

the period F of this function in the frequency domain is just equal 
to the inverse of the time delay t in the time domain. 


This concept can be extended to consider the impulse response 
of the ice layer. For example, in Figure 14. is shown the impulse 
response of a layer 100 cm thick. This function in the time domain 
can be written as 

f(t) = A^6(t) + A^6(t-x) + A26(t-2x) + A3 (t-3x) , 


(A-4) 

where the A's are the amplitudes of the impulses and x is the 
spacing between successive impulses. In Figure 14.,xis approx- 
imately 12 nanoseconds. The impulse response of the ice layer is 
a transient containing theoretically infinitely many impulses of 
decreasing amplitude and alternating sign. This infinitely long 
train of impulses has been approximated by the first four in 
Equation (A-4). A Fourier transform of f(t) then is 
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[f(t)] = F(ju,) = + 

(A-5) 

and with same triganometric manipulation it can be shown that 

|F(jo)) 1^ = + A-j^ + + !\^ + 2 (AqA^ + A^A 2 + A2A3 ) cosojt + 

2(A A. + A,Ao) C0Sai2x + 2(A A«) COSwSt . 

0 ^ 13 ' 0 3 

(A-6) 

If the amplitude values shown in Figure 14. are read off the plot 
and substituted into the above equation, the square of the mag- 
nitude of the transform becomes 

iF(jo)) 1^ = 0.46 + 0.21 coscoT - 0.04 coso)2t + 0.01 coso)3t . 

(A-7) 


But this is just an approximate description of the power reflection 
coefficient for the 100 cm thickness shown in Figure 8. In this 
particular example, the time delay t, in Figure 14., is approx- 
imately 12 nanoseconds, whose inverse, equal to .08 GHz, is the 
period of oscillation of the power reflection coefficient in Figure 
8. corresponding to the 100 cm ice layer. 

If iF(jca) I is subtracted from unity and multiplied by 273, 
the physical temperature of ice in degrees Kelvin, then an approx- 
imate expression for the apparent brightness temperature is pro- 
duced to the extent that the galactic and atmospheric contributions 
may be neglected, i.e. 

T[j (1. - iF(ju) 1^) 273 . (A-8) 

T. 147. - 57. cosojT + 11. cosoi32t - 3. coso)3t 
b 

(A-9) 

This is a description of the apparent brightness temperature curve 
of Figure 9.. for frequencies where the galactic contribution has 
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become negligible. The 100 cm thickness curve corresponds to 
T = 12 nanoseconds, the 10 cm thickness curve to t = 1.2 nano- 
seconds. Furthermore, the time delay x is just the two-way 
propagation time through the ice layer, and therefore directly 
proportional to the ice layer thickness. 
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